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PREFACE 


The idea for a regional seminar on heat transfer and related topics is 
credited to Professor J, Edward Sunderland, North Carolina State University, 
who had both the initiative and courage to organize the first meeting of the 
Southeastern Seminar on Thermal Sciences in April 1965 at Georgia Tech, 

The stated objectives of the conference were; "to stimulate research, 
provide a formal means for presenting research, help research workers in the 
Southeast to become better acquainted with work, colleagues, and resources in 
the geographical area, " Professor Sunderland felt that such a seminar with a 
limited number of formal lectures by prominent engineers and scientists, 
coupled with brief progress reports on current research by people in the south- 
eastern area, would be of interest to the participants. The National Aeronautics 
and Space Administration agreed and fostered the first conference with 
Professor Sunderland, 

Dr. Herbert W, Hoffman served as the general chairman of the Second 
Southeastern Seminar held in Oak Ridge, Tennessee, under the sponsorship of 
the Oak Ridge National Laboratory. Dr. Hoffman reinforced the firm foundation 
established by Professor Sunderland and enlarged the seminar to include topics 
in mass transfer, thermodynamics, fluid mechanics, and energy conversion 
in addition to heat transfer. It was at this second meeting that the official 
name was selected. The adoption of the title "Southeastern Seminar on Thermal 
Sciences" indicated the participants' enthusiasm for the idea originated by 
Professor Sunderland, 

The Third Southeastern Seminar was organized by a committee 
consisting of the organizing committee for the second seminar, and Mr. H. L. 
Atkins of the Marshall Space Flight Center. Recognizing the significance of 
this undertaking, MSFC volunteered to sponsor the third seminar in Huntsville, 
May 8-9, 1967. The organizing committee decided to maintain the format of 
the first two seminars and to present a proposed set of by-laws. These by-laws 
were accepted by the participants, again indicating the continued interest in 
the seminar. 

This volume contains the full text of the invited lecturers and abstracts 
of the progress reports. The group discussions have not been included since 
they were informal. 

The orgamzing committee expresses its appreciation to those who 
participated m the 1967 seminar and to all those who helped in the preparation 
of this report. 
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Definition Of Symbols 


a ,b = constants in general linear viscoelastic model 
mm 

0>j0> t = Jaumann derivative in equation ( 29) 

= nonlinear operator defined in equation (28) 

€ 

G = shear modulus 
g = gravitational acceleration 
h = height of fluid surface 

L = length of tube 
m,n = parameters in power law 
p = pressure 
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Q 

R 

r 

t 

t 

c 

u. 

1 

V 

x,y,z 

Z(a) 

a. 


volume flow rate 
real part of 

radius of tube; outer radius of annulus 

radial coordinate 

time 

characteristic time of viscoelastic fluid 
component of displacement 
velocity vector 
Cartesian coordinates 
Riemann zeta function 

dimensionless rheological parameter in Ellis model and in Spriggs 
model 


/3 

To 

y 


y 




y 
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IJ 


secondary normal stress coefficient 
parameter m several rheological models 
velocity gradient 
elongational velocity gradient 
components of velocity gradient tensor 
complex velocity gradient amplitude 

parameter in viscoelastic models indicating deviations from the 
Weissentaerg hypothesis 

displacement gradient tensor 
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Q = primary normal stress coefficient 


7} = non-Newtonian viscosity 

Tjo = zero -shear viscosity 

7]^ = cross viscosity 

7}^^ = complex viscosity - ~ irj" 

ti'yT]" - real and imaginary parts of the complex viscosity 

7} = elongational viscosity 

K = radius ratio of an annulus 

= time constants in Jeffreys and Oldroyd models 

Aj^ - time constant in Maxwell model 

= component of velocity 

fjL = Newtonian viscosity 

p = fluid density 

T = shear stress tensor 

= components of shear stress tensor 

T = shear stress at the wall of a circular tube 
R 

tj /2 = Ellis model parameter 

= complex shear stress amplitude 
S2 = angular velocity 
0 ) - frequency 

0 )^^ = vorticity tensor components 

t = transpose of a tensor or dyadic product 
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Introduction 


The engineer is familiar with the momentum flux tensor 

T = -)n(W+ (Wj ^) + (2/3/1 - K) (V-V)5 

which describes the mechanical behavior of fluids which are '"Newtonian." This 
expression is inappropriate for polymers, polymer solutions, and considerable 
effort has been expended in the past two decades to develop suitable substitutes 
for the above expression for these structurally complex materials. To date 
not enough has been done to provide a complete fluid dynamical description of 
polymers, but some of the recent results can now be put into a form suitable 
for industrial applications. The purpose of this presentation is to give a 
status report, emphasizing those facets of the subject which seem to have 
potential use by engineers Additional material on theory and applications can 
be found in recent publications [ i, 2,3, 4] . 

We begin by discussing three experiments^ that represent three classes 
of mechanical testing which have been performed to elucidate pol 5 rmeric dynami- 
cal behavior. We see how each of these experiments has led to models capable 
of describing the fluid behavior in situations similar to the original experiment 
Then we show how these experiments have gradually led to the evolution of 
three methods of attacking polymer-flow problems: ( 1) an empirical approach, 
based on die use of constants obtained from viscosity curves combined with 
some experimental data; (2) a model-building approach in which one attempts 
to generate models capable of describing a wide variety of dynamical responses, 
and ( 3) a model -building approach in which molecular -theory results are 
embedded in order to minimize the number of constants needed. 

Throughout we shall restrict ourselves to incompressible fluids and 
isothermal conditions, so that the equations of change are: 

Continuity: 0 = (V-v) (1) 

Motion: ^ ^ 

_ - [V.pw] -Vp -[V.T]+ p'g (2) 


i. Bird, R. B., Lightfoot, E. N., and Stewart, W, E., Selected Topics in 
Transport Phenomena ( 1965) . 
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We shall have occasion also to refer to the Newtonian fluid, for which the 
momentum flux is given here, appropriately simplified for the incompressible 
fluid 


T= _^[Vv+ (Vv)^] = -^xy (3) 

which serves to define the y tensor as the symmetrized combination of velocity 
gradients. We shall also have occasion to refer to the incompressible Hookean 
solid, for which [5] 

T = - G[Vir+ (VT) = -Ge (4) 2 

in which u is the displacement vector and e is the symmetrical combination of 
displacement gradients. Equations (i-4) summarize the starting equations for 
classical fluid dynamics and classical elasticity for materials which undergo 
no volume change on deformation. 

Experiment I: Steady Laminar Tube Flow 


The first experiment we consider is an example of a viscometric 
experiment in which one considers shear stresses in steady flow. 


R. 



We imagine a horizontal tube (Fig. 1) of length L and internal radius 



One measures the volume rate of 
flow Q resulting from an impressed 
pressure difference Ap. For 
Newtonian fluids (N) a linear relation- 
ship is found experimentally between 
Q and Ap, but for polymer solutions 
(P) a distinctly nonlinear relation is 
found ( Fig. 2) . 


FIGURE 1. STEADY LAMINAE TUBE 
FLOW, EXPERIMENTAL ARRANGE- 
MENT 


At low values of Ap the relation is 
linear .and the fluid behaves as though 


it had a large viscosity. As Ap is 

increased, the volume rate of flow increases faster than would be anticipated for 


2, valid only for small deformations 
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a Newtonian fluid' — that is, the "viscosity” appears to be decreasing. Finally, 
at very high Ap there is some evidence that a linear relationship is again 
approached, but with a substantially lower viscosity; it should be remarked that 
this limiting linear behavior is extremely difficult to attain and that at such high 
flow rates viscous heating effects may obscure the phenomenon. An important 
point is that if one plots log Q versus log Ap, there is for many materials a 
large range (sometimes several decades) over which the plot is linear, with 
the lines for various materials having different slopes [6,7] . 



FIGURE 2. STEADY LAMINAR TUBE 
FLOW, EXPERIMENTAL RESULTS 


It is quite clear that in this simple 
experiment Newtonian fluids and 
polymeric materials behave quite 
differently. We now ask how the 
Newtonian behavior can be described 
and whether or not some simple model 
can be evolved to describe the polymer 
behavior . 

Newtonian Fluids. For any fluid, 

equation ( 2) may be integrated to give 

r = Apr/2L, which shows that the 
rz 

shear stress varies linearly with r. 

For Newtonian fluids one inserts into 

this expression the relation 

r = ~ju(dt^ /dr) , which is just 
J^Z z 


equation (3) written appropriately for the geometry under consideration, 
resulting differential equation for 


The 

V (r) can he integrated to obtain the usual 
z 

parabolic velocity distribution, and this in turn may be integrated over the tube 
cross section to give the well-known result of Hagen and Poiseuille. 


Q = 

^ 8pL 


(5) 


which has been adequately verified experimentally 

Polymeric Materials. For polymers one has to try to guess a relation 

to take the place of r - -jj,{dv /dr) which when substituted into 

IlZ z 

T = Apr/2L will yield a velocity profile that in turn will give a reasonable 
Q versus Ap curve. Of the many empirical relationships suggested, a few 
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that have been proposed for polymer solutions are given in "rable I along with 
the corresponding Q versus Ap relationships for circular -tube flow. 

To date, the most widely used model has been the power law, which 
describes the commonly observed lihear log Q versus log Ap plots. Because 
of the ability of the model to describe this behavior and because of the simplicity 
of the model, it has been widely used, often quite successfully. Nonetheless, 
the model suffers from a number of defects which are worth noting. ( 1) Since 
n is usually less than unity, the power law predicts infinite viscosity at 
vanishingly small velocity gradients, ( 2) the parameter m has awkward 
dimensions which depend on n, (3) it is not possible to form a "characteristic 
time" with the model constants — a point that will be discussed m detail later — 
and (4) it has been demonstrated experimentally that if one determines m and 
n from a tube-flow experiment, one cannot use the same m and n for other 
geometries, such as an annulus [i] , a packed column [8], etc. 

Because of the linear log Q versus log Ap plots, it is natural to try 
to make minor improvements on the power law in order to preserve its good 
features and to remove its faults . Three such models have been listed in 
Table I — the truncated power law^, the Ellis model [6], and the Sisko model 
[9] . These three models represent different ways of combining the power law 
and Newton's law of viscosity. We also list the Eyring [ 10] model, inasmuch as 
this expression has been derived from a simple molecular picture. Numerous 
other models have been suggested, but we do not list them here [6,8, 11, 12, 13] . 
Probably more information is available on the Ellis model than for any other 
three-constant model [ 14] . Some sample values of the Ellis model parameters 
are given m Table II [ 15] . 

The models in Table I have been written specifically for a flow situation 

with a velocity profile v = v (r). Eor complicated flow patterns one can 

z z 

generalize these expressions by writing them all in the "generalized Newtoman 
model" form [ 16] : 

' -”^ 1 , 


analogously to equation (3) , in which ?] is a function of the 


T 

ij 


or the 7... 
ij 


1. Spriggs, T. W. , Private Communications 
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But since 77 is a scalar, it must depend on either of these two tensors in such 
a way that only certain combinations of tensor components appear which are 
scalars Now we know that a scalar can be constructed from a vector v by 
forming the dot product of the vector with itself to get (v- v) = 1 ^, which is a 
scalar. 

TABLE I. SEVERAL EMPIRICAL NON-NEWTONIAN VISCOSITY EXPRESSIONS 
AND THEIR INTEGRATED FORMS FOR CIRCULAR-TUBE FLOW. 

= ApR/2L = Wall Shear Stress] 

it 


Model 

Shettr-Btress expression 

Volume rate of 0QW 

Characteristic 

time 

Power law 
(m,n) 

II 

1-1 . 

dv 

z 

dr 

^ (l/n) + 7 ym 

t 

None 

Truncated power law 
< 1 ) 0 . Vo, n) 

dv 

dv 

1 z 

T a -jJo — _ 

rz Yo dr 

... -) 
^ y-dfh«) 

ft) (y 

where r, = i]»Ryo/ 

"r 

L 

$!Sko model 
(Hi. Yfl. n) 

\z~ -•>“ (- 
\\here*c 

e 

1 z 

V /dr < 
z 

n-l\ , 

\ dv 

)-dT 

Yo 

H''] 

t\heic V >.i) and IS Riven b\ ' 

1 1 • ifl) 

1 

Fll IS niudel 
T„ a) 

dv , 

. J- 

tlr no 


T 

F2 

"")v 

•I'lii 




Iiun« (T(|, t,,) 

dv . T 

z 1 1 rz 

-T- - — ‘^inh . 

^ to To 

fr^\\ 

•o Wjl 

[(^) 

slnh (-^ J - co'ih f— ■* 1 1 / 

V'^o/ V'*’9/ Jl 



Similarly, one way to construct a scalar from the tensor y is to form the double- 
dot product of the tensor with itself to get (y:y) = S 13 ■ ?, which is a scalar. 

1 ]ru 

Hence one possible way to generalize models such as those in Table I is by 
writing, for example, 


Power law: 

V = m\ Ji iyiy) \ (7) 
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E1L.S model: 


i 

V 


1 


Vo 


1 + 


(t : t) 



2 


Q!-i 


( 8 ) 


We thus ensure that rj will be a scalar and that the models will behave properly 
under coordinate transformation. As only limited experimental tests of this 
mode of generalization have been made, however, it should be used with caution. 


In conclusion, experiments such as laminar tube flow have led to the 
generalized Newtonian model in equation (6) and to the defimtion of a non- 
Newtonian viscosity r}, which is a function of the shear stress or the shear 
rate. Actually, the tube-flow experiment has today generally been replaced by 
the cone -and -plate viscometer, which gives versus directly, thereby 

removing the need for rather tedious data analysis. Finally, it should be 
emphasized that these non-Newtonian effects can be quite large, with the 
viscosity decreasing by a factor of 10^ to 10® as the shear stress is increased. 


TABLE n. ELLIS MODEL PARAMETERS FOR SEVERAL AQUEOUS 
POLYMER SOLUTIONS AS FUNCTIONS OF TEMPERATURE, [15] 


Solution 

Temperature, 

°C. 

Vo 

poise 

■^1 

^ -2 
dyne cm. 

a 

Range of y 
sec. 

Deviation, 

% 

Max. 

Avg. 



2.2 

49.3 

2.073 

100-1700 

1.6 


1.5% Poly- 

30 

0.58 

148 

2, 165 

100-8000 

7.8 

HI 

ethylene 

40 

0.40 

165 

2.268 

200-3000 

5.1 

1.0 

oxide 

50 

0.235 

214 

2.306 

'200-6000 

6.8 

1.7 


60 

0 ISO 

200 

2.206 

400-7000 

7.3 

1.9 


20 

198.5 

1140 

4.454 

10-900 

8.5 


2 . 0% Hy- 

30 

142.0 

936 

4.068 

10-1000 

10.5 


droxyethyl- 

40 

104,0 

846 

3,920 

10-1500 

10.3 

3.6 

cellulose 

50 

71.0 

760 

3.662 

10-1000 

10.6 

4.8 


60 

39.3 

833 

3.699 

20, 1000- 

6.8 



20 

81.1 

30.3 

2.438 

10-3000 

10.2 


0.5% Hy- 

30 

50.0 

38.0 

2.442 

15-5000 

8.7 


drojqrethyl- 

40 

33.0 

32.8 

2,350 

15-5000 

9.7 

2.0 

cellulose 

50 

21.5 

41.6 

2.486 

30-1000 

5.8 

2.1 


60 

12.3 

47.5 

2.408 

30-2000 

5.8 

1.4 
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Experiment 1 1; Shape Of Free Surface In 
Tangential Annular Flow 


The next experiment we consider is observing the shape of the surface of 
a liquid contained in the annular space between two cylinders, of which the inner 
is rotating with a constant velocity ( Fig. 3) . The Newtonian fluid tends to pile 
up at the outer (stationary) wall because of the centrifugal force. The poly- 
meric materials, on the other hand, do just the reverse — they tend to move 
toward the inner (rotating) shaft and to climb up it; this is called the Weissen- 
berg effect [5, 17] . Here again is a striking difference between the behavior 
of Newtonian and polymeric materials. As in the previous experiment, we show 
how the Newtonian behavior can be explained; then we try to understand what 
modifications in the equations are suggested by the strange behavior of the 
polymeric substances. 
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FIGURE 3 . TANGENTIAL ANNULAR 
FLOW 


Newtonian Fluids. In the system 

under consideration the following 

postulates seem reasonable: v = 

U u 

V -0, V =0, p=p(r, z). Then the 

equations of motion for the Newtonian 
fluid [equations (1,2,3)] become in 
curvilinear coordinates 


? _ 8p 

r “0^ 




(9) 

d 

1 

d 

-] 

(10) 

^d"r 

r 

dr 

(r-j) 

8p 

■ 8z 

Pg 



(11) 


These equations can be solved for the 
velocity and pressure distributions. 

The shape of the free surface is just 
the locus of all points r, z for which 
p(r, z,)= po, the atmospheric pressure. 
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This shape is given by 


h(R) -h(r) = ^ 


I 1 - k2 



( 12 ) 


in which h is the hei^t of the surface and S 2 is the angular velocity of rotation 
of the rotatii^ shaft. 


Polymeric Materials . The Weissenberg effect for polymeric materials 

cannot be explained by the generalized Newtonian models described earlier. 

It seems to be generally agreed now that the origin of the effect is in the 

"normal stresses" , t. t . These stresses are zero for the Newtonian 

tt 69 zz 

fluid and they do not appear anywhere m the equations of motion above. If we 
now postulate that for some reason these normal stresses do not vanish for 
structurally complicated fluids, then the equations of motion contain a few 
extra terms and the result of the fmal expression for the shape of the fluid 
surface which is only approximate is 


h(R) -h(r) - 


1 

g 


R V 


I — 

J r 


dr 

pg 



dr - 




[T (R) -T (r)] 
'■ ,rr ' rr 


(13) 


The dotted underlined terms, which did not appear in equation ( 12) , do mdeed 
contribute to the shape of the surface. 


This analysis suggests that equation ( 6 ) should be generalized in 

some way so as to cause t , and t to be nonzero. About fifteen years 

rr 66 zz 

ago the following suggestion was made almost simultaneously by several 
investigators: 


T.. = 
ij 



1 

2 




i 

3 





(14) 


which is the Reiner-Rivlin-Prager equation [18, 19,20] ; here rj is the non- 

Newtonian viscosity and 17 is the cross viscosity, and both are shear- 

0 

dependent. (Actually they functions of the invariants of the y -tensor. ) 
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When one uses this equation to compute the normal stresses in equation (13) , 

he finds that r and furthermore that r (R) -t (r) is positive if 

rr ee rr rr 

Tj is positive. This is in agreement with the fact that the fluid tends to climh 
c 

up the rotatmg shaft. 

This apparent success of the model gave considerable impetus to its 
further study. It appeared, however, that in rotating-annulus systems of this 
type, ^ but that perhaps (the Weissenberg hypothesis) or 

at least approximately so [21] . 

It IS possible to modify the Reiner -Rivlin-Prager model in order to 
avoid the above objection. One very tentative suggestion is the Bird-Williams- 
Spriggs equation [22]:^ 




+ 0(^y /^t) 

^ IJ 


(15) 


Here is the Jaumann derivative [23] defined in equation ( 29) . It 

represents a time derivative in a coordinate system moving with the fluid and 
rotating with it. In equation ( 15) there are two shear -dependent normal stress 
functions 6 and /?; setting /3= 0 corresponds to the Weissenberg hypothesis 
that in simple shear flow the normal stresses perpendicular to the direction 
of flow are equal. 

The experiment described above was one of the first normal stress 
experiments performed and illustrates strikii^ly the need for formulating 
rheological equations to describe normal stresses properly. It is not, however, 
a suitable experiment for measuring 9. The current methods of precise 
measurement of these normal stress phenomena have been described elsewhere 
[5, 22, 24], 2 


1. Apparently suggested earlier by Ericksen, J. L. , in Viscoelasticity, 
Academic Press, 1960, (see eq. 2.12) . This gives the same results for 
visco metric flow as the simple fluid model (see 2) . 

2. Coleman, B., Markovitz, H, , and Noll, W. , Viscometric Flows of Non- 
Newtonian Fluids, Springer, 1960 
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Experiment III: Transient 
Response in a Cup- Bob System 


We now consider the system in Figure 4. A fluid is placed in a cup, in 
the middle of which a solid cylindrical bob is suspended by a torsion wire. The 
block diagram ( Fig. 4) shows the experimental arrangement with input 
oscillation to cup. The curve (Fig. 4) shows the amplitude ratio for Newtonian 
fluids (N) and polymeric materials (P). One causes the cup to oscillate with 

an amplitude A ( about 1 degree) and a 
c 



CU 03 


FIGURE 4 CUP -AND -BOB SYSTEM 


frequency w . After the initial, 
transients have been damped out, one 
can measure the amplitude of oscillation 
of the bob A. and the phase difference 


(j> between the input and output 
oscillations, both as functions of the 


frequency. Here we consider just the 
behavior of the amplitude ratio • 


For Newtonian fluid Aj^/A^ is 

zero at w = o and w = <» and goes 
through a maximum at the natural 
frequency of the torsion pendulum. 


For polymeric materials, on the 
other hand, the maximum — greater 

than unity — occurs at frequencies higher than the natural frequency of the 
torsion pendulum. Here again we have an example of a striking difference in 
the dynamical behavior of Newtonian and polymeric materials. Let us now 
see how these two kinds of behavior can be described anal34;ically . 


Newtonian Fluids . The response curve for the Newtonian fluid [25] 
can be derived from the set of coupled equations describing the motion of the 
pendulum and the motion of the fluid . We do not go through the details here 
but only state the final result for the special case that the clearance b between 
the torsion pendulum and the eup is small: 
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where M = (v/h^) J l/k and A = 2irIl®Lpb/l, Here v is the kinematic 

viscosity, I the moment of inertia of the bob, k the torsional spring constant 
of the wire, R the radius of the bob, and L its length. 

Polymeric Materials. For polymers [ 25, 26,27] the explanation of the 
large amplitude ratio is sought in the elastic behavior of the fluid. The simplest 
possible suggestion is to combine Newton's law of viscosity with Hooke's law of 
elasticity. If one multiplies Hooke's law in equation (4) hy ij/g and differen- 
tiates with respect to time, then adds the result to Newton*s law of viscosity, 

he obtains (usii^ = y. ) 

ID 1] 

9t,. 

VV-ST'-Wu (17) 


in which and rjo = 2fi, This simple superposition is called the 

Maxwell model [5] . Further differentiation and combining gives 


8r. 


T. + \ 
1] 


ij _ 


8t 


= -Vo 





( 18 ) 


which is called the Jeffreys model. Clearly the process can be continued, and 
one obtains finally 



a 

m 





(19) 


which is called the general linear viscoelastic model. 

Since equation (3) is limited to very small displacements, equation (4) 
must be similarly restricted. If we confine our interest to sinusoidal 
disturbances — as in the experiment under consideration — ^we can postulate that 
the shear stress and the velocity gradient will be undergoing sinusoidal os- 
cillations: 

= e^^^} (20) 

ij ij 
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(21) 


y = m* {7.? 
ij 13 


iwL 
e } 


in which t.® and 7 ^ are complex amplitudes. Substitution of these 
expressions into equation ( 19) gives 

t.?'=-77^^7.® (22) 

IT 13 


where 


17 * = 17 ' - itj" = ;?o 


1 + (xctj) + b2(i^)^ + ... 
1+ ai(iu)) + a^lio))'" + 77. 


(23) 


is the complex viscosity. Clearly the real and imaginary parts 7}* and 77 " 
are functions of co. The quantity 77 ' contains information about the viscous 
response, whereas 77 ” contains information about the elastic response. 

One can now go back to the equations of motion and insert equation ( 19) 

[ or one of its special cases such as equation ( 18) or ( 17)] . These equations 
can be solved and do indeed predict the kind of behavior which polymeric 
systems indicate in Figure 4 [25] . It is further possible to analyze the 
amplitude ratio and phase-shift data to get a tabulation of 77 ’ and 77 ” for a 
polymeric material. A wide variety of experimental arrangements have 
been used by the polymer physicists, and these have been described extensively 
in Ferry's treatise [28] . 


Summary Of The Experiments — The Material Functions 


We have seen how the above-described three experiments [29, 30, 31, 32] 
have led to the definition of various measurable material functions. The tube- 
flow experiment suggests the definition of a non -Newtonian viscosity 77 ( 7 ) . 

The study of the Weissenberg effects has led to the definition of a primary 
normal stress coefficient 0 ( 7 ) and a secondary normal stress coefficient 
j3(y) . The oscillatory flow experiments have suggested the defimtion of the 
complex viscosity 77 =*', with its real and imaginary contributions 77 ' (w) and 
7f"( oj) . A further experiment, which we have not described here, deals with 
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the study of the normal stresses in oscillatory small amplitude motion; this 
leads to the definition of still more material functions [32] 

Since all the experiments we have discussed have involved rather 
complicated geometrical arrai^ements and curvilinear coordinates, it is 
perhaps advisable to summarize the definition of these material functions at 
this point in terms of flow between two infinite flat plates. Such a summary is 
given in Table m [31, 33] . 

TABLE m. SUMMABY OF DEFINITION OF MATERIAL FUNCTIONS 


l. 'steady shear flow i 

2. Small amplitude sinusoidal oscillations 

T= dv /dv 

r: — . — 

X 


. . . tart. 

S V ^ 

T -m' 

\-x ' 


** T ^{d i T ^ 

11 1 11 ^ 

b NoniKil stiess futictioi s 

T - T -O')’ 

XX \y 

n 

Complex viscosity 

L,-! » ri'-ir}" 1 x.y.z 

yx 

T - T — 

yv zz 

b 

Complex normal stress coefficients 

T 0 - T “ 0---e’-i0" 

XX yy 



T « - T “ P- = JS' -1/? 


c. 

Normal stress displacement functions 



^ (d^ -dz} = - 

m. {d -(!} = -/ 1/1=" 
y z 


N 


Note For Newtonian fluids t) = n, n' = M. sU others are zero. 


This summary does not exhaust the possible definitions of material 

functions. These are, however, all the material functions needed for any 

viscometric flow — that is, flows in which there is only one non-zero velocity 

component ( in, say, the i direction) which depends on one coordinate oiiLy 

(but not the i coordinate) , Consider, as an example of a nonviscometric 

flow, steady elongational flow in which v = yz, v = - (i-)yr, 0; this is the 

z r p 

kind of flow which might he approximated m fiber spinning. One then defmes an 


1. For superposition of steady and oscillatory motions, see MacDonald, I. F., 
and Bird, R. B., Journal Physical Chemistry, 70, 2068-2069, 1966. 


18 








elongational viscosity by 


P+ T = rjy ( 24 ) 

zz 

where 7] is Icnown to be dependent upon the elongation rate 7 = di^ /dz. For a 

z 

Newtonian fluid = Bfi, sometimes called the Trouton viscosity. 


Thus it is possible to measure all kinds of material functions. Many 
measurements of 77 are available for many materials, since all one needs is a 
viscometer and since this property has been recognized as useful for predicting 
flow in pipe lines, etc. Frequently ?)* has been measured, since this measure- 
ment is a basic one used for some time in the analysis of polymer structure. 

The primary normal stress coefficient 0 has been measured in a few labora- 
tories, but eq,uipment for trustworthy measurements is limited. Onlj'' several 
isolated attempts have been made to measure the secondary normal stress 
coefficient j 3 , very careful experiments have to be run, and the data have to be 

d d 

analyzed with caution [ 35 ] . To date no experimental curves of Q^, ( 3^,0 , jS 
have been presented in the literature. Also, surprisingly, it seems to be well 
nigh impossible to find trustworthy data on rj in the technical literature. 

What would be highly desir^le are experimental data on quite a few of 
these material functions for a number of polymeric systems. Such data could 
be useful in testing and improving the methods discussed in the next three 
sections . 


Method I: Use Of Non -Newtonian Viscosity Parameters 


In order to describe the mechanical behavior of polymeric materials we 
need some relation for t to use in place of equation ( 3 ) for Newtonian fluids 
or in place of equation ( 4 ) for Hookean solids [ 36 ] . We do not know what this 
relation is, but perhaps we can make some statements concerning what it should 
contain: 


a. Experimental data on 77 indicate that in the limit of zero velocity 
gradient 77 tends toward a limiting value 770. Data on 77' also indicate that 
as the frequency tends to zero, 77' approaches the same limiting value. Hence 
it seems reasonable to suggest that the expression for r ought to contain a 
zero-shear viscosity 770. 
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b. Experimental observations on polymeric materials indicate that 
relaxation phenomena are iniportant. If one suddenly stops the motion of a 
viscoelastic material, there is a relaxation of the stresses in the system and 
this relaxation process can be characterized by a time. We note also that the 
simplest viscoelastic model — that due to Maxwell in equation ( 17) — contains 
a characteristic time. In addition, the Eyring rate theory of viscosity gives 

an expression for tj which contains a characteristic time (Table I) , Hence the 
expression for t ought to contain some kind of characteristic time t^, 

c. Experimental observations on log Q versus log Ap for polymer 
flow in tubes give straight lines over a wide range of variables. The slopes 
of these lines are dimensionless quantities which cannot be formed from a 
viscosity and a time. Hence one might conclude that t ought to contain at 
least one additional dimensionless parameter a. 

The above arguments lead to the conclusion that 


T = VO’ t 


(25) 


where ? is a function that we do not know. All we suspect is that it ought to 
contain certain parameters. 

Now if the ftinction ? were known, then we could presumably calculate 
the various material functions with it to get 


V = ¥y; Vo, t , oi) 

0 = 0(7; %, t , 0!) 


(26) 


If we were to fit, say, data on ?]('y) with a curve containii^ a zero-shear 
viscosity, a characteristic time, and some other parameter a, we should then 
have a set of three parameters to describe the fluid. These parameters, 
although obtained from the measurement of just one material function, should 
nonetheless contaiu information which is also contained in the ? equation. - 
Onee these parameters have been obtained, they can be used in several ways: 
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1. They can be used in the function rj(y) to descibe viscometric flows 

2. They can be used for dimensional analysis of complex flow systems 
[12,15,38,39,40] 

3. The characteristic time can be used to delineate various regions of 
behavior, by comparing it with the characteristic times of the macroscopic 
system [41] . 

4. The zero shear viscosity can be used for getting asymptotic solutions 
for flow problems analytically to tie in with the correlations in step 2 [ 12] . 

5. Attempts might be made to correlate rjo, t , and a with concentra- 

o 

tion, molecular weight, chain length, etc., in an effort to interconnect the 
mechanical and physicochemical properties. 

We see thus that the empirical tj ( y) fimctions discussed earlier may in 
fact have considerably wider utility than might at first be expected — at least 
those three parameter functions which have a characteristic time and a lower 
limiting viscosity. 

A rather complete summary of results obtained by this method is given 
elsewhere [36] . We just mention the experiments which have been performed 
and the models which were used in the analysis: 

a. Correlation of drag coefficients for flow around spheres — Ellis 
model [40] . 

b. Correlation of friction factors for turbulent flow in tubes — meter 
model, Ellis model [38] . 

c. Interrelation of tube flow and annular flow data—Ellis model [ 1, 37] 

d. Correlation of wall corrections for falling sphere viscometer — Ellis 
model [ 15] . 

e. Correlation of friction factors for packed bed flows — Ellis model [39] . 


1. Ashare, E., and R. B. Bird, to be published. 
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f. Correlation of pressure drops through converging sections — Sutterby 
model [ 12] . 

g. Simulation of unsteady state flow in a manometer — Sisko model [41] . 

The last experiment demonstrated that it is probably possible to describe some 
imsteady state motions using only information on provided that the 

characteristic time of the fluid is considerably less than the characteristic 
time of the system (in this case, the time for one manometer oscillation) . 


Method 1 1: Use Of Viscoelastic Models 


A second approach to polymer fluid dynamics [31,42] is to try to guess 
the unknown fiinction in equation (25) . Trying to guess a complicated functional 
relation which will give the proper limiting behavior for all the experimentally 
known material functions is no easy task. We will make no attempt to reproduce 
all the arguments which have led to the kinds of functions which are now 
regarded as useful viscoelastic functions. The discussion here is influenced 
primarily by the trail -biazmg papers of Oldroyd [43,44] . We discuss a special 
case of one of his models — a four-constant Oldroyd model which has consid- 
erable instructional value, since it displays all of the effects we have discussed. 

Let us begin with equation { 18) — the Jeffreys linear viscoelastic 

model — which is known to describe certain low-frequency oscillatory behavior 

of some polymers fairly well. In order to describe nonlinear behavior, we 

replace the linear time operator 9/6t by a nonlinear operator [F : 

s 


r,, + 3^ T. = -rjo (y, + ha ^ T, J 


1] 


e 11 


1] 


e 1] 


(27) 


in which 1> (A 9 /Xd> ~, and the operator 3- is given by reference [33] 

9 £ 
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^ r. = ■— T.. - ( i + €) 
e ij 1] 


TT T., 7i + S T , y, .) 
2 k ik k] jk ki 


t(2 Z 

\m n 


T y l 6. 
mn rnn 1 y 


( 28 )^ 


Here the Jaumann derivative (P3) is 


^ij Bt ^ij^ ^k \ 9x. 2 ^^k^ik\j ^k ^ik^k]^ 


(29) 


The tensor w = (dv /dK - 9r./8x.) is the vortioify tensor and is a measure of 
1] ] 1 1 J 

the local rate of rotation of the fluid. The Jaumann derivative, ^/^t, de- 
scribes the change in as the observer floats along with the fluid and rotates 

with it. Hence, to get equation (27) , one first replaces 9/9t by t in 
the Jeffreys model in order that the latter describe the viscoelastic behavior 
locally, with the rotation and translation separated off. Then one adds some 
nonlinear terms reminiscent of those in the Reiner-Rivlin-Pr^er model — such 
as 2]y y and St. y Finally, two other nonlinear terms are included 

XiC KIJ liC iCJ 

(i-SSr y 6 and ~ SSy y 6.) in order to make t, traceless. 

3 mn nm ij 3 mn nm ij ^ i] 

This results in equation ( 27) , which has been constructed so as to incorporate 
deviations from the Weissenberg hypothesis (that /3 = 0) by inclusion of the 
small parameter e. 


1. The models discussed earlier assume that the stress tensor t is sym- 
metric . One may allow for deviations from S 5 unmetry by introducir^ an 
additional term on the right side of equation ( 28) of the form 




S S r y )d. ] 
m n mn nm ij-" 


in which ^ is an additional constant which indicates the extent to which t 
is nonsymmetrical. Introduction of this extra term contributes to deviations 
from the Weissenberg hypothesis; fiirthermore, it produces deviations from 
the "analogies" [see equations (31,32, 34, 35) and the concomitant discus- 
sion] between the steady shear-flow and oscillatory shear-flow properties. 
To date no experiments have been performed to determine whether or not 
r IS symmetrical. 
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To show how equation (27) is used to obtain explicit expressions for 
the material functions in Table in, we consider first steady shear flow and 
calculate tj, ©, and (3; then we consider oscillatory motion and get rj' and 
7 ]”. 


Steady Shear Flow With v (y) , v = 0, v = 0; 7 = dv /dy. For this 
X y z K 

situation it is desired to compute the six components of from equation (27) . 
These are found to be 


T 

XX 


T + 

yy 


T + 
ZZ 


T 

xy 






^ 2 + 2€\, /2 + 2e\ , 2 


3 3 


mhT 




(30) 


T =0 T = 0 

xz yz 


These equations may be solved simultaneously to give the expressions for the 
individual components of the t tensor. Then from the defimtions in Table III 
one can find the following expressions for the material functions rj, ©, and jS: 


JL ^ 1+ (Ag/?^) (c?47)^ 
770 1 + 


/_2_\ [(1 - (Ag/A^)] (cAj^y) ^ £7 

\ 2 y Tjo 1 + ( cXj7) 2 ey ^7o 


in which 


c = 



IS called the shift factor for reasons which 


will be apparent later. These results show how the various material fiinetions 
depend upon the zero-shear viscosity, the characteristic time Aj, and two 
dimensionless parameters Ag/jq and e, when g = 0, the function /3 is zero 
Note that tins model predicts that the secondary normal stress coefficient is 
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just a constant multiplied by the primaiy normal stress coefficient; there are 
probably not enough reliable data to verify or contradict this at the present 
time [45] 

Oscillatory Small-Amplitude Motion With V (y,t), v = 0, v =0. For 
X y z 

small-amplitude oscillatory motion, equation (27) reverts to the Jeffreys model 

in equation ( 18) . When the postulates in equations (20, 21) are made for the 

Jeffreys model, one obtains a simplified form of equation (23) namely, 


rjQ 1 + iAjO) 


(33) 


When 7]' - 17]" is broken up into real and imaginary parts, one gets 


JT ^ 1+ (V?4) (Aioi)^ 
Vo 

77" _ [1 - (V?4)f (?iiai) 
7]o 1 V (A^OJ) ^ “ 


(34) 

(35) 


The expressions for 6", and © have been derived in detail elsewhere 
[32] for the special case that e =-0. 

The general shapes of the functions in equations (31, 32, 34, 35) are 
correct, but it is not possible to get quantitative fits of the experimental data. 

It IS possible to make order-of -magnitude calculations, however, and these have 
proved usefiil . In addition, it appears from equatiojis (31, 34) that rj/rjo versus 
C 7 ought to be the same as 77’/7jo versus w. Experimentally analogies of this 
type have been found and experimental values of the shift factor c are of the 
order of magnitude one expects from this model. Experimental values of c 
have ranged from about 0.25 to 0. 7 [21, 24, 38,46,47,48] , and the model gives 

c<JW or about 0. 8. Furthermore, comparison of equations (32, 35) shows 

a similar analogy between (c/2) © 7/770 and 77 "/t7o. Analogies of this type 
have also been reported, but they do not seem to be quite so good as the first 
one mentioned. In any case we see here the possibility of interrelating various 
material functions by means of a model. 
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Method 1 1 1: Incorporation Of Molecular Theory Results Into 

Viscoelastic Models 


The iimited success of the Oldroyd models suggested that better results 
might be obtained by working with the nonlinear analogs of viscoelastic models 
of order higher than that of the Jeffreys model. It is possible, for example, 
to generalize equation (19) in exactly the same way that we generalized 
equation ( 18) in the preceding section. The model that results has too many 
constants to be of any value. Hence the House molecular theory of linear 
viscoelasticity [49] was used to evaluate the parameters and b^ in 

terms of a single time constant A. In that way the two-constant Spriggs-Bird 
model was obtained [50]. 

In an effort to improve on the latter, Spriggs [33] reformulated and 
extended the two -constant equation just mentioned, and it is his results that we 
choose to summarize here, since they tend to tie together many of the ideas in 
^ this chapter. The starting point for the discussion is equation ( 19) , which we 
wish to generalize to handle nonlinear problems. First equation ( 19) can be 
written in an alternate form as a sum of Maxwell models: 


CO 

T = J t(p) 

P=1 



(36) 


(37) 
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SPRING-DASHPOT REPRESENTATION 



and the constants a and b can be related to the time constants \ used 

mm p 

here. If one now puts in ' X = one gets a good approximation to the re- 

P 

suits of the Rouse molecular theory 

Spriggs suggested that the Rouse model thus formulated be extended in the 

following way: { 1) replace 9/9t everywhere by (equation 28) , (2) set 

cx. ^ 

A. = X/p . In this way one obtains a model containing four constants:^ 


rjo = zero-shear viscosity 

A = characteristic time (Rouse's longest relaxation time if a = 2) 
e = a dimensionless parameter which is small and accounts for 
deviations from the Weisseriberg hypothesis 
O' = a dimensionless parameter which governs the slope of the log rj 
versus log y line in the "power-law" region, usually between 2 
and 4. 


This model has been extensively studied by Spriggs who has derived expressions 
for the material functions and given limiting forms which are computationally 
easy to use. It seems to give material functions for rj, 17', and 17" which are 
capable of fitting the experimental data quite well; the 0 data do not seem to be 
described as well as the other properties. 


Tlie model also predicts exactly the same analogies as are predicted by 
the four -constant Oldroyd model with exactly the same shift factor c. Hence 
all the material functions ( 77 , rj\ r}'\ 0, |S) can be obtained from Spriggs's 
results for rj* and 77": 


v'(^) _ V ^ 

1 77"(w)_Z(2a) 

' p=l p (p + (Aw)2) 

Aw T7o Z( a) 

(Aco) 2 y X 


Z( 0 !) 2a. '2a 

' ' p=i p (p + (Aw)2) 



(38) 


(39) 


in which Z( a) = ^ n is the Riemann zeta function. 
n=i 


1. Another nonlinear extension of the general linear viscoelastic model was 
made independently by Roscoe [58] , he did not, however, incorporate any 
molecular theory results into his model. 
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Spriggs has also shown how to relate his four-constant model to the 
Ellis model discussed, by noting that both models exhibit linear log ri versus 
log y behavior beyond about /Vcy = 1 . For large y the asymptotic forms of 
the two models are 


Ellis: T 

yx 


Spriggs: 


i-a 



y a 


i-a 1 


no 


IT 


, , /a+ 1 \ 

Z(a)2a sin^— ^ ttJ 


X (cA) y 


(40) 

(41) 


From this it is evident that a has the same meaning in both models and 
fiirthermore that 


cX = 



7T 


Z(a)2a sin 



a-1 


(42) 


This gives a means of estimating A from Ellis model parameters, or vice 

versa. For a = 2 and e = 0, the time constants are related by 

A =s 0.56(?7o/Ti) • Hence the characteristic time of the Ellis fluid is of the same 
2 

order of magnitude as the largest relaxation time of the Rouse theory. Since 
some work has been done on the temperature dependence of the Ellis parameters 
[ 15] , this information can be transferred to the Spriggs parameters hy means 
of equation (42) . 


Summary Of Analytical Methods 


The preceding methods must be regarded as tentative, since as more 
data accumulate and various new models are devised and tested, undoubtedly 
better procedures can be established. In any case considerably more is 
known about viscoelastic models today than was known twenty years ago, and 
the three methods summarized here show some promise. 

The kind of rheological models that are to be used in any given problem 
depends on the variety of phenomena vMch must be Interrelated, on the visco- 
metrie data available for the fluid, and on the degree of complexity of the flow 
configuration. 
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For polymer structure studies, small amplitude oscillatory responses 
are measured; this kind of "mechanical spectroscopy" is well described by 
Ferry [28] . For the interpretation and design of those experiments, the 
general linear viscoelastic model is quite adequate. Furthermore, considerable 
progress has been made in understanding these problems from a molecular 
point of view [5i, 52] . 

For engineering problems involving steady flow in closed channels, the 
simple generalized-Newtonlan models are probably sufficient. Even some 
unsteady problems can probably be worked, with only viscometric data, 
provided that the time constants for the macroscopic system are larger than 
the characteristic time of the fluid. For heat transfer problems in steady state 
flow, probably viscometric rj data will be sufficient. It will, however, be 
required to have information concerning the temperature dependence of the 
viseometric parameters; this information is known for the power-law model and 
the Ellis model for a few fluids [ 15] . Most of the existing analytic solutions for 
heat transfer problems have just assumed that the rheological parameters are 
insensitive to temperature [53, 54, 55,56, 57] , A notable exception is Turian’s 
recent study of viscous heating in a cone-and-plate viscometer [ 15] , 

For systems mvolving viscoelastic recoil phenomena, large amplitude 
oscillatory flows, spinning operations, and normal stress phenomena, nonlinear 
viscoelastic models are required [3] , It would be helpful to be able to describe 
these highly nonlinear phenomena by equations which contain only three or four 
constants at the most and to be able to determine these constants in one or two 
well-defined experiments . For that reason, the discussion here has been 
restricted entirely to models containing small numbers of constants rather than 
functions. Such models show more promise for use by the engineer, (For a 
comparative summary of rheological models see; Spriggs, T. W. , Huppler, 

J, D,, and Bird, R, B., Trans. Soc, Rheol,, 10:1, 191-213, 1966, This 
paper includes both "differential models" and integral models. ") 

Illustrative Examples 


Laminar Flow of an Ellis Fluid in a Thin Slit 


Problem: Derive an expression for the flow rate of an Ellis fluid in a 
sHt of width W, thickness 2B, and length L. Let the z axis go down the middle 
of the channel so that the side walls oi the slit are located at x = iE and 

y = ±1 W. 
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Solution: The equation of motion for any fluid in a slit is 
, dr 

f. ^ ^ 

dz dx 

We may integrate this using the boundary condition that t =0 at x = 0; 

xz 

r = (Ap/L)x 
xz 

Introduction of the Ellis model into this result gives 


dv 
z 

dx 


Jl 

^0 


+ 

I Lt. 


cn-l 


fApx 


which can be integrated with the boundary condition that v = 0 'at x = ±B to 

. z 

give 




= 4* nu- + 


Br^/^o \2 r ^ ■ OH-1 


where g = x/b and n = ApB/LTjt. Then integration over the flow cross section 
gives ® 

Q= 2WB f V d^ 

or 

Q 

for flow in a plane slit. 

We may adopt this result to flow in an annulus formed by two tubes of 
radii kR and R, provided that k > 0,6. We are then dealing with a slit for 
which 2B- ( i - x)R and W= 27rR[l - i (1 - k)J . We can in addition multiply 

the final result by the fimction by which the analogous 

Newtonian slit formula is "corrected for curvature" in order to describe annular 
flow. The final result is 


2WB'^Ti 
a_ 

Vq 


n + 


CW-2 


n 


O' 
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Q = 


7rR3(l-K)2[l-l(l-K)][l+ 1 (l-/c)2+ 

_ 



Tl 

3 


X 



n + 


1 

ctt-2 



where n= (1 - f£)RAp/2Lri. 

3 

Laminar Flow of a Polymer Solution in an Annulus. 

Problem: For an aqueous 3.5% CMC-70-medium-solution Ellis 
parameters were determined from Fredrickson’s [ 1] tube flow data to be 

•> 7 q =22.7 poise ^ 

Ti = 1, 520 dyne cm.” 

2 

O' = 3.0 

Find the flow rate (expressed as 4 Q/ ttR^) of this fluid' in an annulus with radius 
ratio K = 0,624, when ApR/2L= 5 x iO^ dyne cm. 

Solution: Use the result of illustrative Example 1 mth 1 - k = 0.376, 
i - ( 1 - k)= 0.812, and 


n= (0.376X5 X 10V( 1,520) = 1,24, 


Then 


4Q _ 4(0, 376)2(0. 812)(1. 002) (1,520) 
'W (22.7) 


1 (1.24) + i- (1.24)3 


= 24.3 sec. 


This result agrees very closely with the experimental value obtained by 
Fredrickson [ 1] , 
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stretching of a Viscoelastic Fiber. 


Problem: A polymeric material [34] is being drawn out as a fiber. A 

relation of the form of equation (34) is desired in order to relate the tension, 

p + T , to the rate of elongation of the fiber, y = dv /dz; that is, we wish to 
zz z 

know the elongational viscosity. Further, we want to relate this quantity to the 

Shear viscosity (non-Newtonian viscosity) , which is normally measured with 

a viscometer, such as a cone-and- plate instrument. 


Solution: We consider a coordinate system oriented so that the central 
axis of the fiber is the z axis, and we assume that the flow is steady. 


elongational flow with p = yz, vr = — — yr, y0= 0. This assumed velocity 

Z /j 

distribution satisfies the equation of continuity exactly. It is further assumed 
that the stretching is sufficiently slow that the term [V« pw] is not important 
in the equation of motion; also external forces are neglected. These 
assumptions lead to the fact that the shear stresses are independent of position, 
a point we shall use subsequently. The first thing we have to do is to summarize 
ihe input information for equation (27, 28, 29) for the assumed flow pattern. 
Because the problem involves cylindrical coordinates (r, 6, z) we have to know 

various operations specifically Vv and {vVr} in cylindrical coordinates. 

These have been tabulated by Williams and Bird [32] , and we use those results 
here. We find that 


Vv = (Vv) 


t 



7 


0 0 


y = Vv + (Vv) 


t 


-1 0 

0 ~1 

0 0 


L 


+ i 

0 

0 

+2 


7 


0 ) = Vv - ( Vv) = 0 


rr 




0 

0 


69 


0 

0 


; , {vVt} = 0 


zz 
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{r.y} = 







“ 

-T 0 

rr 

0 


1 

0 

0 


0 

11 

• 

• *N 

i 

0 

i 

0 

o 

o 

1 

+2t 

zz 


0 

-r 

0 

4 

_ 


When this information is put into equation ( 27} , the three nonzero components 
of that equation become 


T -A^(l+e) [-r 7-4-(-t - 2t )y] 

rr ^ '■ rr 3 rr 60 zz 


=-ho {-Y - ^ (i-*- e) [r^ - 4 (67^)]} 


[-Tg^r -T ^-'^rr - V ^^zz^^ ^ 
=-^7o {-7 -A^d+e) [y^-~ (67^)]} 


— I _ 

T - Aj (1 + e) [+2 t 7 - — t - r „ + 2 t ) 7 ] 
zz ^ zz 3 rr 00 zz 

=-Voi+^y - ?i2 ( 1 + e) [47^ - (67^)]} 


Solvir^ for t and t gives 
zz rr 


2r = 
rr 



1 - Aa( 1 -i- e)7 
1 - Ai( 1 + e)7 


At the outer edge of the fiber there should be no normal forces, so that 

p + T =0, whence -r =+ _1 t . Hence 
rr ’ ^ rr 2 zz 
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P == I ^zz = 


i - Asd-s- e)y 
1 - Ai(l+ e)7 

Comparison of this with equation (24) gives 


V = 3^0 


1 - ^2(1+ g)y 
1 - Ai(l+ €)y 


Since 74 > 7^, this predicts that tj will increase with the rate of elongation. 
Note that r) decreases with y for the same viscoelastic model. For 
Newtonian fluids 74 = 74 =0, and we get the classical result that rj = Bjx, 

We can now try to interrelate rj and rj by using this result and that in 
equation (31). We define two functions S and T thus 


S(y) 



and 


T(y) = - 


dy 



These functions can be tabulated when non-Newtonian viscosity data are 
available. 'These functions can also be calculated from equation (31) . Then if 
we set e = 0 and eliminate 74 and 74 among the expressions for ’rj’, S, and 
T, we get 


JL 

3r]o 


1 



S( 1 - S) - T 1 /i 

S -T J V 2 


S -T 
T 


1 - 



S - T 
T 


which gives an estimate of tj from viscosity data. There do not seem to be 
data available to test this relation. 

Tube Flow of an Oldroyd Fluid. 

Problem; Obtain the volume rate of flow for the model given in n 
equation ( 27) (the four- constant Oldroyd fluid) [ 42] , 
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Solution: As pointed out earlier, the equation of motion gives for tube 

flow 

Prom equation (31) we find that 


^ _ 1+ n(o?4y)^ 

i+ (c?47)2 ^ 


where n = Combining these we get 


r= 1+ n(c\^Y)^ 

‘^P 1 + ( cTy.j^'y) 2 

This could be integrated again to get the velocity profile, but we want only 
Q versus Ap. Hence we write the volume rate of flow Q in this way by two 
integrations by parts: 



where y - -dp /dr| _-p;7-p is related to the pressure drop Ap by 

Jtl- Z -Tv 

HcT 

t)o \2h/ \i+\l s 

where When the integral in Q is evaluated, we get 

22i£-=ri. 1 f„(n-x„)l He: 

where 

(n;Xj^) = 1 n 3 x ^2 _ 3 ^ 2 ^^, _ 3^,^^ - 1 ) ( 2 n - 1 ) 


xin(i + Xj^> -Ixj^ ^rn4)* f«" 
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From the above equations for Ap and Q the quantity has to be eliminated 

to get the Q versus Ap relation. This can he done numerically, and 
dimensionless plots are available both for circular tubes and for plane slits. 

The results have been compared with experimental data for earboxymethyl- 
celluiose, sodium polymethacrylate, and cholesteryl butyrate [ 42] . 

Relative Magnitude of Normal and Shear Stresses. 

Problem: In simple shear flow what is the maximum value of 
]t^ |/|t^ I for the Oldroyd model in equation (27) when e = 0 ? This gives 

some idea as to the relative magnihides of the normal and shear forces. 


Solution: Let R be the desired ratio. Then from the definitions in 
Table HI and the results in equations (31, 32) 


R = 


fyyi 

It 1 



(^ 

1 + 


Setting dR/d 7 


0 gives 




\ , and so 


R 

max 


(Aj - A2) 



As AgAi goes from 1 to 1/9, goes from zero to Js/3, 
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AXIAL ANNULAR FLOW 

-^AA/^- A/AAAA- A/VW 



PRESSURE-SENSING DEVICES 
REAP SAME AT A AND B 


DIE STREAK 


UH 



FLUID FLOWS 
WHEN PU«P IS 
TURNED OH 



PRESSURE-SENSING DEVICES 
READ DIFFERENTLY AT A AMD I 


SECONDARY FLOWS 
-A'VWVW' 

ROTATING DISK 
BEAKER 





[see c.t, hill, 

J.D. HUi»PlER, 
AND R.6. eiRIK 
CHEH. ENGR.SCI., 
( 1966 )] 


COMPARATIVE RESULTS OE NEWTONIAN AND 
POLYMETRIC FLUIDS 
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FOUR TYPES OF MOLECULAR THEORIES 


1. SUSPENDED PARTICLES EINSTEIN 

(ANNALEN DER PHYSIK 1906, 1911) 

2. PEARL NECKLACES ZIMM 

{JOURNAL OF CHEMICAL PHYSICS, 1956) 

3. ACTIVATED STATE EYRING 

{JOURNAL OF CHEMICAL PHYSICS, 1936) 

4. NETWORK LODGE 

(TRANS. FARADAY SOCIETY, 1956) 



BEAD-SPRING MODEL 

FLEXIBLE JOINTS 

SPRINGS ARE "GAUSSIAN 
CHAINS" 

SPRING CONST. = ^ 

MEAN SQUARE END- 
TO-END DISTANCE 


I 


ZIMM 
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JUNCTIONS DON’T FLUCTUATE; HAVE FINITE LIFETI/VIES 
CHAINS ARE G^\^SSM; FORM COHERENT NETWORK 
AFFINE DEFORMATIONS 
NO POLYMER- SOLVENT 
INTERACTION 

USE RUBBER THEORY RESULTS 


LODGE 




ELLIS MODEL (3 CONSTANTS) 


J, ± 

V 


1 + 



a-n 


1. ZERO-SHEAR VISCOSITY 77 

'0 

2. CHARACTERISTIC TIME r) It,, ^ 

'0 1/2 

-3. CONSTS. HAVE SIMPLE DIMENSIONS 
4. SIMPLE, YET USEFUL 


SPRIGGS- BIRD MODEL ( a-2, €'0) 
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POLY- 


[-CH2-CH2-] 
[- CH2 - CHCl -] 
[-CH2-CH-] 

6 

[- Si{CHj)2 -0-‘ 


ETHYLENE 

VINYL-CHLORIDE 

STYRENE 

DIMETHYLSILOXANE 



STEADY LAMINAR TUBE FLOW 



MEASURE : VOLUME RATE OF FLOW 
Q, AS A FUNCTION OF PRESSURE 
DROP AP 
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TANGENTIAL ANNULAR FLOW 




OSCILLATOR VISCOMETER 



MEASURE: AMPLITUDE 

'^ATlOAgog/A^^jpASA 

FUNCTION OF THE 
FREQUENCY. 


-BOB ALSO OSCILLATES 
WITH FREQUENCY w, BUT 
WITH AMPLITUDE Ar 


DIFFERENT FROM A 


CUP IS MADE TO OSCILLATE WITH FREQ.w 


'BOB 

CUP. 




OSCILLATORY VISCOMETER 



I 


NATURAL FREQUENCY OF BOB 


DISK IN CYLINDER EXPERIMENT 
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7 EXPERIMENTS TO TEST MODELS CONTAINING 
■»7 AND A CHARACTERISTIC TIME 


1. TURBULENT TUBE FLOW 

2. LAMINAR ANNULAR FLOW (AXIAL) 

3. FLOW IN A CONVERGING CHANNEL 

4. FLOW IN A PACKED COLUMN 

5. DRAG COEFF I C I ENTS OF S PHERES 

6. SPHERES FALLING IN CYLINDERS 

7. OSCILLATING MANOMETERS 
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(FREDRICKSON) 

(McEACHERN) 

(SUnERBY) 

(SADOWSKI) 

(SLATTERY) 

(TURIAN) 

(BIERY) 
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FRICTION FACTOR, f 


METER'S FRICTION FACTOR 
CORRELATION 
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LOW TEMPERATURE THERMODYNAMICS 


By 

Richard E . Sonntag 
University of Michigan 
Ann Arbor, Michigan 


Introduction 



Thermal analysis of processes involving simple substances and their 
mixtures at near-ambient conditions has not commonly presented any unusual 
problems, because in this range the behavior of such systems can be predicted 
or correlated by simple and familiar models. In recent years, however, it 
has become necessary to consider the thermodynamic behavior of these sub- 
stances at extreme conditions of both temperature and pressure, each extreme 
having its own class of problems that does not lend itself to as simple an 
analysis . 

The purpose of this discussion is to consider the low temperature 
region, in which problems of thermod3mamic phase behavior become of funda- 
mental interest and importance. Two classes of processes in which phase 
equilibrium problems arise will be discussed, with particular emphasis on 
the low temperature system hydrogen-helium. 


Condensation Problem 


To illustrate this problem, let us begin by considering the simple 
process of cooling a pure gas (for example, nitrogen or oj^gen) in a heat 
exchanger. The thermodynamic properties and various phase boundaries of 
most pure substances in the cryogenic region are reasonably well known at 
low and moderate densities, and for a few substances at high density. Thus, 
the gas behavior during cooling and the conditions at which various phase 
changes will occur can be predicted with a high degree of confidence. Of 
course, heat transfer and flow considerations are vitally important to such an 
analysis, but not to the purpose of this discussion. 


55 



This condensation process becomes much more complex if a mixture 
instead of a pure gas is being cooled, and xs much less capable of being 
satisfactorily analyzed at the present time. Knowledge of die thermodynamic 
behavior of the gas mixture (enthalpy, for example) , requires experimental 
measurements, calorimetric orP-V-T-Y. Except for the low density 
region, these have been almost non existent in the low temperature range 
until the past few years. Consequently, analytical correlations of this type are 
not well developed. The phase equilibria can, in principle, be predicted from 
pure component phase behavior and P - V - T - Y behavior of the individual 
phases, but again, experimental data for the latter are very sparse at present 
and efforts in predictions of this type have not been generally successful. 

Certain special problems can be correlated by analytical techniques, 
however, and help to give an insight into the more general problem. As an 
example of such a simplified phase equilibrium problem, consider the 
selective freezing of component "A" from a binary gas mixture of "A" and "B" 
being cooled to a low temperature. This process may be used in the puri- 
fication of gas "B," or the situation may simply arise in cooling the gas 
to low temperature for liquefaction. In either case, it is important to know 
the amount of "A" being solidified from the mixture onto the heat exchanger 
tubes . 


This situation can be classified as a binary solid vapor phase 
equilibrium problem, and commonly the condensed phase can be assumed to 
be pure "A. " It is not the purpose of this presentation to develop an analysis 
of this problem. Nevertheless, a brief discussion of a set of assumptions 
commonly made m developing a simplified model will help to point out the 
regime in which such a model is or is not reasonable. The requirement for 
the' equilibrium of "A" between the two phase at a given T, P can be expressed 
in terms of the fugacities as 


^OL ^ |GAS 
V A 


( 1 ) 


The solid phase fugacity at T, P is given in terms of the fugacity of saturated 
”A*' at the same temperature T and saturation (sublimation) pressure 
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The fugacity of component "A" in the gas mixture can be expressed in terms of 
the P - V - T - Y behavior of the mixture. Thus, the problem can, in princi- 
ple, be evaluated for the composition of the gas mixture at equilibrium. 

Let us now make a set of three assumptions regarding this phase 
equilibrium: 


1. The right-hand side of equation 2 is negligibly small; i.e. , 

SOL ^ SAT 
^A ^A 

SAT 

2. Pure saturated vapor "A” at T, P^ behaves as an ideal 

- sat „ sat 

gas;i.e. , f^ = P^ 


3. The gas mixture at T, P behaves as an ideal gas mixture; 

i e T = Y P 

0 


Substituting, equation (1) becomes 

P = Y P ' 

A A 


(3) 


which is the equilibrium equation commonly used in analyziiig familiar 
thermodynamic problems of this type (e.g. , saturated air -water vapor 
mixtures at near ambient conditions) . Of the assumptions made in developing 
this simplified model, the third assumption is generally the most serious. 
Careful consideration must be given to the regions in which equation ( 3) is used 
to predict equilibrium composition, because it may be in error by unexpectedly 
large amounts. A convenient parameter for evaluating the system departure 
from ideal behavior is the ratio of the actual mole fraction Y^ to that given 
by equation (3) , Y^ , or 

o 


e = 





SAT 


(4) 


A number of solid-vapor systems have been studied experimentally in 
recent years, and compared with analyses that include the corrections and 
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utilize accurate equations of state to predict the mixture P - V - T behavior. 

An example of this type of investigation is the CO 2 - Ng system [ i] , with COg 
solidiftring from gas mixtures of the two components. The apparatus required 
for experimental measurements of this type is shown schematically in Figure 
5, and a set of results, experimental and analytical, in Figure 6. It is seen that 
the departure from ideal behavior becomes extremely large, with the gas-phase 
concentration of COg being 250 times that predicted by equation (3) at the 
highest pressure studied. 

It is also seen that the various analytical equations used m the more 
general analysis give very good predictions of the composition within the range 
of validity of the equations themselves. The analytical solution is involved and 
tedious, even for a simple phase problem as posed here. For more complex 
equilibria, predictions become even more difficult, and in all cases rely very 
heavily on the accuracy with which the P - V - T - Y behavior of the phases 
can be represented. 


Tank Pressurization Problem 


The storage or transfer of a cryogenic liquid by pressurizing with 
an "inert" gas, commonly helium, leads to another consideration of binary 
phase equilibria m the low-temperature region; this is a liquid-vapor 
equilibrium problem. Under most conditions, there will be only a small 
amount of helium in the liquid, and the vapor will be mostly helium. Knowledge 
of the phase compositions at different temperatures and pressures is a 
necessary prerequisite to analysis of a transfer process. 

The general lack of present ability to predict even this problem has 
led to a number of experimental studies of liquid-vapor equilibria in binary 
systems. The hydrogen-helium system is one of partictilar interest in the 
field of cryogenics, and has been investigated over a fairly extensive range of 
temperature and pressure during the past several years. The majority of this 
work [2] was done using normal hydrogen for convenience, and comparisons 
made using parahydrogen [3] . 

A schematic diagram of the vapor-recirculation apparatus used in this 
research is shown in Figure 7, and a portion of the results in Figures 8-11. 

Of particular interest in this study were the solubilities of helium in the 
liquid phase as a fiinction of T and P, given in Figure 9 . 
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Measurements of liquid and vapor compositions in this system have more 
recently been extended to higher pressures [4] , using the same type of technique. 
The results shown on the plot of Figure 12 have a shape in the vicinily of 28®K 
that is characteristic of a system exhibiting a high - pressure phenomenon 
commonly termed gas-gas equilibria. This occurs where the system critical 
locus changes to a positive slope on a P - T projection, and is believed to exist 
for the hydrogen-helium system, although not directly confirmed. 

Another phenomenon studied in this research was the phase-densily 
inversion locus, that condition under which the vapor (helium rich) phase 
becomes more dense than the liquid (hydrogen rich) phase, and settles to the 
bottom of the containing vessel. This locus is shown in Figure 13. 

In connection with the problem of liquid hydrogen storage, it has heen 
proposed to store cryogenic hydrogen in a slush form, mixtures of the solid 
and liquid phases. The technique appears promising at present, and has two 
significant advantages over liquid storage: higher density, with resulting 
decrease in storage volume; and decreased storage losses, resulting from the 
heat of fiision of the solid and heat capacity of the subcooled liquid. 

Consequently, a study to evaluate the lower -temperature liquid -vapor 
and solid-liquid-vapor equilibria in the hydrogen-helium system was under- 
taken, and has recently been completed [5] . A glass dew -bubble point 
apparatus, Figure 14, was used to determine the various phase boundaries 
in this region. Results of particular interest are the liquid-phase compositions 
in the liquid-vapor region. Figure 15, and the liquid and vapor compositions 
on the three-phase surface, shown in Figure 16. Research on this type of 
system is being extended at the present time to evaluate the effects of ambient 
and other heating on the storage of slush mixtures. 
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SOME ASPECTS OF AERODYNAMlfTlEA'TlNG 


By 

R. W. Truitt 

North Carolina State University 
Raleigh, North Carolina 


Introduction NTO»1g0gg 


The problem of continuum aerodynamic heating on a body involves the 
solution of the boundary layer over the body. Generally speaking, the 
solution of the boundary layer over a body depends on first, determining the 
local inviscid (or aerodynamic) flow over the body. The simplest body is, of 
course, that of an infinitely thin flat plate at zero angle of attack and, as a 
first approximation, it becomes a trivial aerodynamic problem since one may 
assume with good approximation that the pressure gradient is zero in all 
directions. On the other hand, for certain combinations of Mach numbers and 
Re5molds numbers, the boundary layer thickness on the flat plate may cause 
the flow to ’’feel” a body of finite thickness such that the assumption of zero 
pressure gradient is no longer valid. The solution to the "flat plate" boundary 
layer then becomes considerably more complex and the classical boundary- 
layer-shock-wave interaction phenomenon must he considered as a dommant 
factor in the aerodynamic heating problem. If the boundary layer flow over a 
simple or complex body shape, either two-dimensional or axisymmetric, is to 
be calculated, the local inviscid flow field must be found. In the cases in which 
aerodynamic heating is involved, this means that one is confronted with the 
rather formidable task of determining the local inviscid flow field where the 
local pressure, density, temperature, etc., are varying from point to point 
over the body. Unfortunately, the flow and thermophysical properties of the 
fluid are strong functions of the altitude, speed of flight, the kind of gas, and 
the attitude and shape of the body. Numerical techniques of varying complexity 
and accuracy are available, not only for solving the inviscid fLow field but for 
the local boundary layer as well. No completely satisfactory solutions for 
either the inviscid flow field or the boundary layer have been found that consider 
the actual non-equilibrium flow over a typical blunt-nosed fli^t vehicle at 
hypersonic speeds. However, some real progress has been made that allows the 
approximate solutions for some important limiting cases. 
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It is the purpose of this paper to review the factors involved in the 
calculation of the aerodjmamic heating over a typical blunt-nosed body at 
hypersonic speeds, with particular attention to the problem of determining 
the local mviscid flow field under the most realistic flow condition, i, e , 
non-equilibrium flow. It should be emphasized that this is undoubtedly the most 
important factor in the aerodynamic heating problem, since the local freestream 
conditions of pressure, temperature, density, degree of dissociation and 
ionization must be known before any attempt can be made to solve the 
corresponding boundary layer problem. 

Once the local conditions have been determined, the general governii^ 
equations of the boundary layer must be cast in a form that considers the 
viscous non-equilibrium flow of a reacting gas under gradients in the flow 
direction of pressure, temperature, density, dissociation, ionization, etc., 
and for wall conditions that must involve the chemical character of the wall 
(material, surface condition, catalysis, etc.). At the present time, some 
rather sophisticated numerical solutions for the non -equilibrium hypersonic 
laminar boundary in air have been calculated for flat plates under special wall 
conditions. The results found are valuable, but in no way represent what could 
be classified as satisfactory for practical aerodynamic bodies. No satisfactory 
solutions for separated, transitional, or turbulent boundary layers have been 
found for non-equilibrium hypersonic flows . 


Perhaps the classical case that embodies all of the various conditions 
of flow regimes and thermophysical conditions is that of the flow over a blunt- 
nosed body at hypersonic speeds ( Fig. 17) . Here one sees the characteristic 
bow shock wave standing off from the body a certain distance, d. The 
formation of the shock represents the mechanism of heating the gas as it flows 
from the relatively cold freestream, experiencing an abrupt rise in pressure, 
temperature and density, a decrease in velocity and a decrease m Mach 
number. Generally speaking, the abrupt rise in pressure, density, temperature, 
etc. IS accompanied by abrupt changes In the character of the gas. That Is, 
before the gas passes through the shock and indeed as it passes through the 
shock and experiences the abrupt rise in pressure, temperature, etc. , in one 
to two mean free paths, the gas is said to have certain "active'’ and "inert" 
degrees of freedom. For a monatomic gas, the active degrees of freedom are 
those in translation, i. e. there are 3 active degrees of freedom, and the ratio 


of specific heats 




where z is the number of (classic) 


degrees of freedom. Thus for a monatomic gas under, say, room temperature, 
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7 = 1+ — = 1.667+. For a diatomic gas at room temperature, there are 3 

O 

degrees of freedom in translation, and 2 in rotation that are "active" degrees 

2 

of freedom and thus y = 1+ = 1.4. Now the inert degrees of freedom for 

5 


the diatomic gas are those in vibration, dissociation, ionization, etc., in which 
energy is taken from the gas in order to cause these "excitations." For a 
monatomic gas there is no vibration or dissociation mode, but the gas may 
undergo excitation of the planetary electrons and experience high excited states, 
single ionization, double ionization, etc. Counting each excitation as a degree 
of freedom, it is clear that 7 i for very high temperatures. 


When the cold gas passes through the shock, the active degrees of 
freedom are those in translation (for a diatomic gas, the 2 rotational degrees 
of freedom) . The high temperatures generated then begin to activate the 
"inert" degrees of freedom. Now, whereas the active degrees of freedom are 
adjusted in one or two, or a few collisions (mean free paths) , i. e. , are 
adjusted "instantaneously, " the inert degrees of freedom require several or 
many collisions to come to so-called equilibrium conditions. Since the gas is 
moving and time is Involved in bringing the gas to equilibrium in the vibrational, 
dissociation or ionization mode, it is clear that there is a distance over which 
the gas is in a non-equilibrium state. The extent of the non-equilibrium regime 
downstream of the shock is one of the classical problems in the hypersonic 
flow over a blunt body and the nature and extent of the non-equilibrium flow 
determines, to a large extent, the character of the mviscid flow over the body, 
i. e. , the local inviscid flow, which in turn determines the condition and nature 
of the local boundary layer. 


The inviscid flow over the body may be thought of as the flow through 
a certain convergent-divergent nozzle (or streamtube) from a fictitious 
reservoir represented by conditions in the stagnation region of the blunt body. 
Clearly, then, the character of the flow, 1 . e. , frozen, equilibrium or 
non-equilibrium, along the convergent-divergent streamtube, is a sole 
fanction of tiie initial conditions at the stagnation point. The initial conditions 
at the stagnation point are, of course, a hinction of the altitude and Mach 
number and, perhaps to some extent, the shape of the body. At any rate, for 
ture hypersonic speeds one would anticipate the gas, if diatomic, to be either 
partially or fully dissociated and partially or fully dissociated and partly 
ionized; if a monatomic gas, it would be electronically excited and probably 
ionized. For either gas, the rate equations for the various reactions will 
require a certain time (or distance) to reach equilibrium and the problem is to 
determine if the flow has had time (or enou^ distance) to reach equilibrium 
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at the stagnation point and, if it has, or has not, to determine what the condi- 
tions are at the stagnation point, i.e., the degree of dissociation, ionization, 
etc. 


In the following, it will be of interest to examine briefly the various 
real gas problems in the hypersonic flow over a blunt body that are necessary 
to determine the local inviscid flow from which the boundary layer is determined 
and thus the skin fricition and heat transfer, i.e. the aerodynamic heating. 

The first to be considered is the non-equilibrium flow downstream of the normal 
portion of the bow shock. This allows the determination of the conditions m 
the "reservoir’' or stagnation region. Then the classical "equilibrium" nozzle 
flow is summarized for a reacting gas and some of its characteristics are dis- 
cussed. A more general analysis of the equilibrium and frozen thermophysical 
properties are discussed in which both dissociation and iomzation are considered. 
The frozen and equilibrium cases are important, since they represent the 
"bounds" of the actual gas state. That is, the actual gas flow is generally in 
the non-equilibrium condition which lies somewhere between the two extreme 
states of frozen (zero reaction rate) and equilibrium (infinite reaction rate) . 

This will be carried out for the simplest possible gas, hydrogen, and then 
results for the so-called "air-like" ideal dissociating gas will be indicated, 
both for nozzle flow and the blunt body problem. The implications of these re- 
sults in the calculation of the boundary layer skin friction and aerodynamic 
heating will be indicated. Finally, a review of the boundary-layer-shock -wave 
interaction problem will be given and its effect on aerodynamic heating will be 
discussed. 


Normal Shock Analysis 


In order to treat the boundary layer over a blunt body (which is shown 
to an exaggerated scale in Figure 17) it is first necessary to determine the 
flow conditions through the bow shock wave. Since the external inviscid flow 
just outside the boundary layer represents the local "freestream" conditions 
needed for the calculation of the local boundary layer, the flow conditions aloi^ 
the stagnation streamline must be determined For given freestream conditions, 
the solution for the flow through the normal portion of the bow shock gives the 
first step in determining the flow properties along the stagnation streamline. 

An example of a normal shock analysis is outlined in T^le IV Details 
of this analysis are given in Reference i. This analysis is carried out for 
argon and considers atom-atom and atom-electron collisions and takes into 
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TABLE IV. NORMAL SHOCK ANALYSIS 



A + A^A + A'*’+e 

A+e^A'^+e+e 

(1) 

Mass: 

pu - p^U 

(2) 

Momentum: 

2 ,,2 
p + pu * P^U 

(3) 

Energy: 

i + u^/2 - U^/2 

(4) 


where the enthalpy, i, of the gas is 

i - (5/2) (k/m)(T^ + e T^J + (k/m) e (5 

p/p = (k/ir) [T3 + e Tg] 

T/- (A/2) { (5/8) [1 - (4/5)^ a~z)/X]^'^ + (3/8) - (c/X) 1 - £ Tg*(7) 
P* - (Po/Pion^/ ^ 1 - (2 /X)It/ + e T^*] - (c/X) ) (8) 

X * mU"/2 k (9) 



T^/T 


ion 


* ^e^^ion 


if 


P 


P/Pfon 


( 10 ) 
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TABLE IV. NORMAL SHOCK ANALYSIS (Concluded) 


u {dc/dx) - (p/m) { (1 - e) [(T^^c/'^a) + 2} exp (-W/Tg) X 

[(1 - e) - ( £ A 3 exp Cl/T^*))/Ta3^2j + 

3/2 

e A2 Tg l<Texc''’^c>+ ('■^exc'"*^eHa - c) - (U) 

( E A3 exp (1/T/n/T^^^^l } 

in which, for argon, = 0.1488 x 10 A 2 = 0 3759 x 10*^^, and 
A 3 = 0 3451 X 10'^® and is the energy of the first excited 
state 

* *3/2 

[(3/2)T^* + 1] { (1 - £)T^ ^"■^exc''\^ t ^ ^ 

- ( c A3 exp (1 /t/))/t/^^^] + c ^ 

exp (-Tejcc/Tg>t<^ - - ( e A 3 exp (l/T *) ) /T^*3/2 j Agc/^ X 

{ ( e A 4 /T/) log^ (A 3 //77 ) + (1 - eHAg + (T^* - 0 016)^] IX 

(Tg* - Tg*) ( 12 ) 

where, for argon, = 0.1004 x 10"^, A^ = 0.7064 x 10”®, A^ = 

0.1600 x 10"^®, A 7 = 0 5015 X 10”^^. Ag = 0 1510 x 10^ 


X ~ A„ 


\r ^ - (^A)]^/[F(T^*,e ) +G(T/,e)]}de 

o 

where Ag = 0.2230 x 10'^\ 


03) 


\ - k T /n^^^ » 5 

D a 


(14) 
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account the difference between the atom and electron temperatures during this 

non- equilibrium flow process. The integration of equation ( 13) (Table IV) 

allows the determination of the flow properties throu^ the shockwave. 

Figure 18 shows the variation of the pressure, p, density, p, degree of 

ionization, e, the atom temperature, T, and the electron temperature, T , 

© 

through the shock in the form of dimensionless parameters. The example 
shown is for a Mach number of approximately 27 and reveals that the shock 
thickness is 0. 19 cm (0.075 in, ) . The shock thickness is defined here as the 
distance required for the flow to reach equilibrium conditions and is of key 
importance in determining the flow conditions at the stagnation point on the body. 
The relation of the shock flow to the blunt body stagnation point can perhaps be 
better understood by referring to the schematic curves in Figure 19 , This is 
essentially the same character of the solution shown in Figure 18 but depicts the 
critical relations between the shock detachment distance, the shock thickness 
(shown as the relaxation distance) and the stagnation point. The equilibrium 
flow region is actually an iscntropic deceleration of the flow from ]ust down- 
stream of the shock to the stagnation point. The conditions at the stagnation 
point are strong flinctions of the freestream conditions, the kind of gas, the 
shock detachment distance and the relaxation distance. It is entirely possible, 
under certain conditions, that the relaxation distance is greater than the shock 
detachment distance. It is obviously important to know when and if this situa- 
tion occurs, since the entire flow from the stagnation point downstream along 
the body depends on the initial "reservoir,” or stagnation point, conditions 

Stagnation Point Conditions 

Once the stagnation conditions are known, then the thermophysical 
properties of the gas can be calculated to get the "starting” conditions of the 
gas state in the stagnation region. Table V gives the governir^ equations for 
deterraing the properties of molecular hydrogen, H 2 , undergoing tiie various 
reaction processes of dissociation, ionization, and electron attachment. 
Molecular hydrogen is chosen here, since it is the simplest diatomic gas that 
still retains the essential physical and chemical features of an air-like gas 
For given thermodynamic cond^ions one can thus determine the mole fraction 
of the various species (ions, H , electrons, e~, atoms, H, molecules, H 2 , 
and atoms with an electron attached, H ) . Other important parameters can 

be found, namely, the equilibrium speed of sound, a , the frozen speed of 

© 

t 

sound, a^, the frozen ratio of specific heats, the equilibrium ratio of 
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TABLE V. GOVERNING EQUATIONS FOR MOLECULAR H2 UNDERGOING 

REACTION PROCESSES 

BASIC EQUATIONS 
Hg— /3[ H* Hg +■ H- 

CONSERVATION OF H NUCLEI (PROTONS)- 

2=)9,■^;83+2^l^ + ^g 

CONSERVATION OF ELECTRICAL CHARGE* 


0 “ A ■ Ai ■ A, 


'8 


H H*+ 


+ «■ 


^ (T)= AM 

^ ' A.AgP 
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A,A„P 


MOLE FRACTION* X 


'8 '"8 
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lA, 
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« RT 


RT 




^*1 .- * [r^ 
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TABLE V. GOVERNING EQUATIONS FOR MOLECULAR H2 UNDERGOING 

reaction PROCESSES ( Continued) 
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TABLE V. GOVERNING EQUATIONS EOR MOLECULAR Hg UNDERGOING 

REACTION PROCESSES (Continued) 


GIBBS-HELMHOLTZ ^ E « -T 


F «-kT ln^"/N!] 


E»Nkf 
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TABLE V. QOVEBNINQ EQUATIONS FOR MOLECULAR Hj UNDERGOING 

REACTION PROCESSES (Concluded) 


DISSOCIATION RANGE, 

8/3,+ 28/% =0 

8J + RT8/>.8m , 8P-n 
T mp m p 

8m + Z(m-m.)8^i = 0 

8^4 ISA 2 8/3b 8P AHm 8T -O 
lA /8, P RT* 

AHl 8/% + E^,C. 8T - SJ 8 p = 0 

Hf ‘ F, P 
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Normal Portion of 



HGUHE 19. THE VARIATION OE TEMPERATURES AND DEGREE OF 
IONIZATION IN THE SHOCK DETACHMENT REGION (not to scale) 



specific heats, y , etc Since dissociation is essentially complete before 

ionization begins, it is possible to treat the governing equations in "steps" 
as for the dissociation range. 

Typical results for hydrogen are shown in Figures 20 through 24 

Figure 20 verifies that the process of dissociation is essentially complete 

before ionization begins Figure 21 shows the variation with temperature of 

the ratio of the equilibrium and frozen values of the speeds of sound, the ratios 

of specific heats and the specific heats at a constant pressure, for a pressure of 

one atmosphere Figure 22 shows the average molecular weight of the gas as 

the temperature is increased from 10^ to 10® ®K, for one atmosphere. At the 

lowest temperature the gas is a diatomic one and as the temperature increases 

the "plateau" represents the attainment of complete dissociation to a monatomic 

gas. Further increase in temperature results in the attainment of complete 

ionization of the hydrogen gas. Figure 23 shows that, for a given pressure, an 

increase of temperature yields osciUatory-type variations in the equilibrium 

ratio of specific heats that finally approaches the classical monatomic value 

of i 667, as the temperature is increased to large values. The effect of 

decreasing the pressure is to "squeeze" the curve. It is interesting to note 

that the envelope of the curves approaches y = 1, at high temperatures 

0 

Figure 24 shows the equilibrium speed of sound as a function of temperature 

and the attainment of complete dissociation and ionization are evident as 

"plateaus” in the curve. Figure 25 shows typical results for y for air and it 

0 

is clear that the character of this curve is quite similar to that found for 
hydrogen (Fig 23). 

The foregoing discussion summarizes the type of analysis required to 
determine the thermophysical properties of a typical air-like (diatomic) gas 
under either equilibrium or frozen conditions in the stagnation region. Although 
tedious, the calculations are quite straightforward but require a detailed 
knowledge of the quantum mechamcal properties of the gas which can be 
obtained from the NBS Tables or other similar sources,. 

Inviscld Blunt-Body Real-Gas Flow 


For simplicity, it is usually assumed that the flow at the stagnation point 
is under equilibrium conditions As the flow accelerates from the stagnation 
point around the blunt body, the flow will probably depart from equilibrium as 
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the gas begins to expand and cool. Whether or not the gas species recombine 
to form atoms and molecules (and thus release heat) , or "freeze" at some 
fixed composition, must be determined by a detailed numerical analysis. 

In order to calculate this flow over the body it is advantageous to know, 
if possible, the pressure distribution along the body surface (stagnation stream- 
line) . Fortunately, the pressure distribution over a blunt-nosed body can be 
determined by a combination of Newtonian Impact Theory, Prandtl-Meyer flow 
and Blast Wave Theory [2] It is particularly important to note that this 
pressure distribution (shown schematically m Figure 26) is essentially 
insensitive to whether the flow is frozen, equilibrium or non-equilibrium. On 
the other hand, the flow properties along the body are very strong functions of 
the flow conditions (frozen, equilibrium or non-equilibrium) as well as the 
stagnation conditions . 

As noted previously, it is helpful to consider the flow along the body 
surface as analogous to the flow through an equivalent convergent -divergent 
nozzle (Fig. 17) . With this conceptual model, it is simple to identify the 
stagnation point on the body as a fictitious reservoir, the sonic point ( SP) on 
the body which corresponds to the "throat" of the equivalent nozzle, and a 
downstream station, A, as a supersonic station of the nozzle. Since the 
pressure distribution is known from hypersonic theory, an isentropic analysis 
in terms of pressure ratio will yield the area ratio, A/ A*, similar to that 
depicted in Figure 27. Using this approximation for the area ratio of the body 
surface streamtube, it now remains to calculate the real -gas flow through 
this fictitious nozzle, the results of which should approximate closely the 
inviscid flow properties along the body surface 

Table VI outlines a nozzle-flow analysis for a triply-iomzed monatomic 
gas under equilibrium stagnation conditions. This analysis is one of the 
simplest examples of a reacting gas flow that still retains the essential 
physical features of a typical reactii^ monatomic gas. Details of this analysis 
are given in Reference 3. Interpreted as a blunt-body problem as described 
above, the stagnation conditions are such that the argon gas is 50 percent 
triply ionized (Fig. 28) . As the flow accelerates from the stagnation pomt 
and expands (to larger values of A/a*) the gas recombines first to a fully 
( 100 percent) doubly ionized gas (first transition region) . As the flow 
continues to expand downstream, the gas continues to recombine to a fully 
(100 percent) singly ionized gas (second transition region) 

Figure 29 gives the complete solution for the flow properties for the 
given stagnation conditions including the variation of the pressure, p, the 
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TABLE VL NOZ ZLE -FLOW ANALYSIS 


governing equations for the Isentroplc, equillhrlum flow of a triply- Ionized 
monatomic gaa in a convergent-divergent nozzle are- 
Lav of Mass Action 

nj.n^/n^_ 3 ^ = 2 ( 2 jm^kT/h^)^^^tQ^(Ay)/Q^(Aj._j^)]exp[-{ ,r-l,2,3 ( 15J 

Conservation of aectrlc Charge 

i rn = n ^ 16 ) 

r-1 ^ e 

Conservation of Nuclei 

i n = n, ( 17 ) 

r-o *■ ^ 

Equation of State 

P » kT(n„ + i n ) (18) 


Conservation of Energy 
= V^/2 + E + P/P 
Specific Entropy 

S “ E/T •* M /p «= constaint 
Specific Internal Energy 


E - (kl^/p) ^/5T 


Conservation of I ’.ass 


pAV = 


where 


<1 = E n -fln((2Hm,kT/h^)^^^ Q (A )/n ] + 1 - 1 / k?) 
r^o ^ A err r J 

♦ n^{ln[2(2inB^kT/h^)^'^^/n^] + ij 


“ So(A^) + g^(Aj.) exp(- j|) 
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density, p, the velocity, v, and the temperature, T, along the body (equivalent 
nozzle) . An interesting property of the flow, which is typical of the flow of 
any reacting gas over a blunt body, is the sudden drop in temperature along the 
body corresponding to the location of the transition regions ( Pig. 28) . Although 
the density is not influenced by these sudden drops of temperature, the pressure 
is influenced considerably as noted by the large pressure drops in these 
transition regions. 

It IS important to note that these pressure drops may have important 
implications in the boundary layer development along the body. The large 
pressure drops may tend to thin out the boundary layer in the transition region 
but the less favorable pressure gradient downstream of this region may create 
a situation that is destabilizing and may cause the boundary layer to tend to 
separate 

The preceding equilibrium solution for a triply-ionized gas flow m 
a nozzle is typic al of that for any reacting gas , in which only one reaction is 
occurring under given thermodynamic conditions. This is best demonstrated 
by examining the results of a nozzle flow analysis for an ideal dissociating 
gas. Figures 30 through 33 represent the solution for the flow of Lighthill's 
ideal dissociating gas (IDG) [4] in a typical hypersonic nozzle. The curves 
denoted with the letter "E" represent the equilibrium solution similar to that 
outlined for the ionizing gas, except that only one reaction, the dissociation, 
is considered for the IDG, whereas three reactions, separated by transition 
regions, had to be considered for the triply -ionized case. In addition to the 
equilibrium solution, the solutions for the hypersomc nozzle flow for the 
frozen (denoted by "F") and the non- equilibrium (denoted by ”N") cases are 
given in Figures 30 through 33. The frozen solution is actually the well- 
known "cold gas" analysis (i.e. , a non-reacting nozzle flow, with the same 
reservoir conditions as the other two cases: T = 5900 ®K, p = i. 15 x 10^ atm.) 
The non-equilibrium solution was first discussed by Bray [5] and uses 
Freeman's [6] rate equation in conjunction with the Lighthill ideal dissociating 
gas. 


An examination of the three solutions given in Figures 30 through 33 
show a large difference in the values of the degree of dissociation, the 

local temperature parameter, T/T^, the local velocity parameter, u/u^, and 

the local Mach number, M, for given values of nozzle station, A/A*. The 
dotted line shown in Figures 30 through 32 is Bray's simplified approximate 
solution to the non-equilibrium case. 
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FIGURE 21. EQUILIBRIUM TO FROZEN FLOW RATIOS VERSUS TEMPERATURE 





FIGURE 22. AVERAGE MOLECULAR WEIGHT VERSUS TEMPERATURE 
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p =.001 atm. 



PIGURE 25. EQUILIBRrUM RATIO OF SPECrPIC HEATS VERSUS TEMPERATURE FOR AIR 
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FIGURE 28 . RESULTS OF NOZZLE FLOW ANALYSIS AS APPLIED TO BLUNT BODY 
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FIGURE 29 FLOW PROPERTIES FOR 



BODY ANALYSIS 




FIGURE 30. a VERSUS A/A^> IN A HYPERSONIC NOZZLE 
o 

(Stagnation conditions- = 5900 °K p = 1. 15 x 10^ atm) 
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FIGURES!. T/T^ VERSUS A/A* IN A HYPERSONIC NOZZLE 


(Stagnation conditions; T =5&00“K p = 1, 15 x 10^ atm) 

0 o ' 

(Solid lines - '*exact'* solution, Dotted line = "approximate" solution) 







FIGURE 32. U/U^ VERSUS p/A* IN A HYPERSONIC NOZ ZLE 
(Stagnation conditions: T = 5900 p = 1. 15 x iO^ atm) 
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FIGURE 34. GEOMETRY OF THE SHOCK WAVE AND THE BODY SURFACE 

(S = Shock Wave) 
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FIGUEE 36. VARIATION OF v/u« ALONG THE AXIA L STREAML] 
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FIGUEE 38. VAEIATIONOF p/p ALONG THE AXIAL STEEAMLINE 
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FIGURE 39. variation ALONG THE BODY SURFACE 
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As might be expected, the non -equilibrium solution, which is essentially 
an exact solution, is the most difficult case requiring a numerical solution. In 
all cases shown, the non-equilibrium solution hes somewhere between the two 
"bounds” (i.e , the frozen and equilibrium cases) . It is important to note 
that the reservoir conditions are assumed to be under equilibrium conditions as 
"starting values" for the three cases treated. Under'these assumptions, it is 
clear that, if one interprets A/A* as the location along the blunt body surface, 
the local temperature, degree of dissociation, velocity, local Mach number, 
etc. , are strong functions of the character of the flow (frozen, equilibrium, 
non-equilibrium) . The fact that these values represent tlie local inviscid 
conditions that would be used to calculate the boundary layer, indicate the 
importance of determining the most realistic inviscid flow conditions Taking 
the non-equilibrium solution as exact, it is apparent that the incorrect assump- 
tion of either frozen or equilibrium flow could introduce large errors in the 
corresponding boundary layer solution. 

If one interprets these results in the design of a hypersonic wind tunnel. 
Figure 33 shows that the conventional "cold gas" analysis (frozen case) would 
predict, say, a Mach number 20 for a given A/A^', whereas the actual flow 
( non-equilihrium with no boundary layer correction) would produce a Mach 
number of only 10. A detailed study of these results will convince even the 
most skeptical of the importance of obtaining the non-equilibrium solution for 
reactii^ flows at hypersonic speeds 

Non-Equilibrium Flow 


As already indicated, the inviscid real gas flow over a blunt body at 
hypersonic speeds is a strong function of the shock detachment distance and 
the starting or "reservoir" conditions (i.e. , whether the flow is frozen, 
equilibrium or non-equilibrium at the stagnation point) . It was shown that 
the attainment of equilibrium flow in the stagnation region is not always 
guaranteed, although the nozzle solutions discussed did assume equilibrium 
conditions in the reservoir (or at the equivalent body stagnation point) . It is 
clear from the foregoing that the non- equilibrium condition in the stagnation 
region malces the flow over the blunt body considerably more difficult to 
calculate than for the frozen or equilibrium cases . The flow over a blunt- 
nosed body was calculated for the frozen, equilibrium and non-equilibrium 
cases but for the same freestream, hypersonic conditions [7] . The model 
gas was essentially the Lighthill IDG where the Freeman rate equation was 
used for the non -equilibrium case. 
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The inverse method was employed wherein the shock shape is assumed 
and the shape of the corresponding blunt body is determined together with the 
solution of the inviscid flow aloi^ the stagnation streamline. The results are 
shown in Figures 34 through 41. Figure 34 shows that, for given freestream 
conditions and given shock shape, the shock detachment distance is largest for 
the frozen flow case and smallest for equilibrium conditions. Figures 35 through 
38 show the variation of the flow parameters along the stagnation streamline 
from the station just downstream of the normal portion of the shock ( = 0, 

where R is the radius of curvature of the blunt body) to the body stagnation 
o 

point Note that in all these figures the frozen flow case terminates at the 
stagnation point corresponding to the shock detachment distance, 

5 = ~ 0- 08, whereas the equibbrium flow solution terminates at 6 ~ 0.063. 

The non -equilibrium solution shows the shock detachment distance to be less 
than that for the equilibrium solution. Figure 35 shows that the degree of 
dissociation at the stagnation point for non-equilibrium flow is less than that 
for the equilibrium case The corresponding values of the velocity parameter, 
v/u^, are plotted in Figure 36 for the three flow conditions. Figure 37 reveals 
rather large difference m the local temperatures, although the equilibrium and 
non-equilibrium stagnation values are nearly identical and are much lower 
than the temperatures for the frozen flow case. The local density ratios plotted 
in Figure 38 show that the equilibrium and non-equilibrium values are nearly 
the same at the stagnation pomt and are considerably larger than the corre- 
sponding values for the frozen flow. 

Starting with these values of the parameters in the stagnation region, 
the local flow aloi^ the blunt body was calculated for the three cases and the 
results plotted in Figures 39 through 4i. Figure 39 shows the variation of the 
velocity ratio, u/u^, for the frozen and equilibrium cases; the non -equilibrium 

solution coincides with the equilibrium solution. The variation of the degree 

of dissociation, a , is shown in the lower curve of Figure 39. The local 
0 

body temperatures are seen (Fig 40) to be the lowest for the non-equilibrium 
solution and highest for frozen flow The reason that the equilibrium solution 
yields higher temperatures farther along the body, as compared to the non- 
equilibrium case, is due to the large number of atoms recombining and con- 
sequently releasing more heat Figure 41 depicts the pressure distribution 
for the three cases The curve denoted as is the Newtonian (impact) 

pressure distribution which is plotted for comparison. The pressure distri- 
butions are quite similar in character with only slight differences in the 
magnitudes of the pressure ratios. This is convincing evidence of the fact. 
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as indicated earlier, that the pressure distribution is relatively insensitive to 
real gas effects. It is clear from this example that the local values of Resmoids 
number, Mach number, temperature, etc., are strong functions of the flow 
Conditions (frozen, equilibrium or non-equilibrium) and would therefore result 
in significantly different values of boundary layer skin friction and heat transfer 
for the three cases. 


Laminar Boundary Layer In Hypersonic Real Gas Flow 


The calculation of the boundary layer over a blunt body in hypersonic 
flow depends on an accurate knowledge of the local inviscid flow properties just 
outside the boundary layer. An. inviscid flow analysis similar to that outlined 
above will yield these required local flow properties. The analysis of the 
boundary layer for the flow of a reacting gas over a blunt body is quite complex 
and, as yet, no satisfactory solution has been found, in the discussion to follow, 
a relatively simple solution of a laminar boundary layer for a reacting gas will 
be described. In particular, the results of an analysis for the frozen and 
equilibrium laminar boundary layer flow of a singly-iomzed monatomic gas over 
a flat plate will be discussed The results of this simple solution should never- 
theless reveal the essential features of the laminar boundary layer of a reacting 
gas flow. 

The governing two-dimensional, laminar boundary layer equations for 
a reacting gas are summarized in Table VII, namely, the equation of state, 
the conservation of species, the continuity equation, the x- and y- momentum 
equations and the energy equation. The details of the derivation of these 
equations are given in Reference 8 . A modified Karman-Tsien-Lees solution 
is found [8] for the frozen and equilibrium flow over a flat plate for both 
an ideal ionizing monatomic gas and an ideal dissociating gas. Typical re- 
sults are shown in Figures 42 through 47 for the ionized monatomic gas. It is 
important to note that the frozen and equilibrium solutions represent the upper 
and lower "bounds" so that the more realistic non-equilibrium solution lies 
somewhere between these two limits. 

Figure 42 gives the velocity profiles for Mach numbers 5 and 10 for 
three local "freestream" values of degree of ionization, namely, 0.6, 

0 933, and 0.997. The Blasius solution (M^ = 0) is plotted for comparison. 

The results presented are for the "cold wall" case smce the degree of ioniza- 
tion at the wall, = 0 . The important features to be noted in Figure 42 are 
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the large growth of the boundary layer with increase in Mach number, the 
influence of the local xnviscid degree of ionization on the velocity profile, and 
the nearly linear character of the velocity profile at Mach number iO. Figure 
43 shows the differences in the temperature profiles for the frozen and 
equilibrium solutions for a moderate hypersonic Mach number. Note that the 
frozen case predicts much higher temperatures near the wall than for the 
equilibrium case . Figure 44 gives a composite plot of the frozen temperature 
profiles for Mach numbers 0, 5, and 10 for different values of freestream 
ionization For a given Mach number the main differences in temperature 
occur in the outer region of the boundary layer. 

One of the most important parameters in hypersonic boundary layer 
flow is the boundary layer displacement thickness. That is, the effective body 
shape IS found by adding the local boundary layer displacement thickness to the 
corresponding geometrical ordinate of the body. Thus, the present flat plate 
has an effective "parabolic" shape that will experience a corresponding hyper- 
sonic pressure distribution. Figure 45 shows the very marked growth of the 
boundary layer with increase in Mach number It is also clear that increasing 
the degree of ionization in the local inviscid flow tends to "thin out" the 
boundary layer. Finally, Figures 46 and 47 show the differences in the 
character of the ionization profiles for frozen and equilibrium flows for Mach 
number 5 for = 0. 60. Of particular interest is the "knee" in the equilibrium 

curve and the absence of one for the frozen case which corresponds to a "hot 
spot" in the boundary layer for the frozen flow ( Fig. 43) . 


Boundary Layer-Shock Wave Interaction 


At hypersonic Mach numbers the character of the laminar boundary 
layer flow over a flat plate'results in high skin friction and heat transfer rates. 
Even higher aerodynamic heating rates occur for a turbulent boundary layer. 
Details of the theory of laminar and turbulent boundary layers and the techniques 
for calculating the aerodynamic heating over bodies are given in Reference 9, 
including a full discussion of the maximum aerodynamic heating condition at the 
body stagnation point Another critical aerodynamic heating situation that may 
rival or exceed the stagnation point heating rate is when a shock wave in 
hypersonic flow interacts with a turbulent boundary layer. This problem has 
been treated for the case of protuberances in Reference 10 and will be 
discussed briefly in the following. 
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Figure 48 shows a plot of the ratio of the local heating rate, h, (due to 
boundary layer-shock interaction) to the flat-plate heating rate, h , (with no 

interaction) for the case of a cylindrical protuberance mounted on a flat plate 
at Mach number 3.51. It is clear that the detached bow wave which "wraps" 
itself around the cylinder interacts with the turbulent boundary layer creating 
a maximum heating region on the plate near the normal portion of the bow 
shock ( Fig. 49) . A side view schematic of the flow situation is shown in Figure 
50 where the shock wave interacts with the turbulent boundary layer (on 
centerline) giving rise to an abrupt thickening of the boundary layer and the 
formation of a lambda-footed shock. This abrupt thickening of the boundary 
layer acts as an equivalent wedge angle, Q , with an attached shock wave and 
creates a pressure jump and a corresponding increase in aerodynamic heating 
in this region, A detailed mathematical analysis of the flow model is given in 
Reference 10 and the results are plotted in Figures 51 and 52. The pressure 
rise predicted by the theory is plotted in Figure 5 1 along with experiment as a 
function of the product of Mach number, and Reynolds number, Rey^. It 

is clear that the agreement between theory and experiment is excellent. 

Figure 52 shows a parametric plot of the maximum heating rate, h , as a 

function of the ne-v^ hypersonic turbulent boundary layer parameter, Xj_. The 

correlation with experiment is considered satisfactory. The Important and 
rather surprising point is that, under these interaction conditions, an extreme 
"hot spot" is generated ahead of the protuberance that can be many times 
greater than that of the simple undisturbed turbulent boundary layer for the 
same flow conditions. For example, Figure 48 shows the maximum local "hot 
spot" to be more than 7 times hotter than the heat transfer rate for the 
undisturbed turbulent boundary layer for a moderate Mach number of 3.5 1. It 
is to be noted that the heating rate for the undisturbed turbulent flat plate 
boundary layer at hypersonic speeds can be of such intensity that the ablation 
of steel can easily occur. For true hypersonic flow the "hot spot" can approach 
values of h/h^ much greater than 10, so it is easy to see that this represents 

one of the most critical design conditions for hypersonic vehicles. 

Flow Separation 


Another severe aerodynamic heating problem occurs when the super- 
sonic flow separates upstream of a geometrical flow deflection and reattaches 
downstream (Fig. 53) . It is clear that the separation region as well as the 
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reattachment region can both be treated by an analysis siroilar to that of 
Reference 10. It is interesting to note the differences between the separated 
flows for laminar and turbulent boundary layers ( Fig. 54) . The pressure 
rise oan be predicted by the model used in Reference 10 by using the skin 
friction and boundary layer displacement thickness expressions for laminar 
and turbulent flows. 

Boundary layer separation on delta wings has been studied [11] . The 
separation phenomena is shown in Figure 55 and can be treated by a method 
similar to that of Reference 10 as regards the formation of the lambda-footed 
shock and the calculation of the pressure rise and heat transfer. Note that a 
vortex has formed downstream of the point of separation. This additional flow 
parameter complicates the flow pattern considerably and, as yet, has not been 
treated analytically. 

Real Gas Effects In Shock Interaction 


The analysis used m Reference 10 can be modified to take into account 
the real gas effects. This can be accomplished by taking into account the real 
gas effects across the oblique shock portion of the lambda-footed shock 
(Fig. 50) . Figures 56 through 59 give the pressure, temperature, Mach 
number and Reynolds number as a function of the angle of incidence, a, which 
corresponds to the ai^ie 6 of Reference 10 (Fig. 50) . These figures show 
clearly the difference between the assumption of a perfect gas (frozen flow) and 
an equilibrium flow. These results are based on the analysis of References il 
and 12. To use these results in the analysis of Reference 10, it is necessary 
to derive an expression for the ratio of the skin friction coefficient due to 
interaction, c^, to the skin friction coefficient without interaction, c^ , as a 

fiinction of the corresponding ratios of temperature, T/T^, Mach number, 

M/M^, and pressure, p/p^> The desired analysis is outlined in Table Vm 

using the notation of Reference 10. Note, in this analysis, that the assumption 
that the velocity ratio, V/V^ — 1, which was made in Reference 10, is not made 

here. The last expression for Table Vin can be used in conjunction 

o 

with Figures 56 through 58 to calculate for the frozen and equilibrium 

o 

cases, provided one knows a, the boundary layer displacement thickness 
(referred to as 6 in Reference 10) . It is clear from the results shown in 
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TABLE Vni. C. AND ANALYSIS 

I f 

0 


c, = C(Rey,r 

'w 


c, = C(V/VJ (yoyyO^KReyJ 


0 ' ' W' 


Cf = C(p/pJ (T„/T„) (V/V„) (ReyJ 

% = C(l/4 )(Rey„) 


C = VL= 1+ Pr°(7-I)/2M‘ 


(v/y,) 


. / \* + md + n) 

VV 'p/Po> (v/v„) 

^ T(l+Pr'’(y-l)/2Mf) , 

^VL>Co= Z JT— 1 - (X/T)(M,/M) 

. ' T(l + Pr (y-l)/2M ) 

1/2 

(V/VJ = (M/M„)(T/TJ 


= <T/TJ 


- [l + m(l + n)J 2 + m- 1 2 [l +m(l+n)] I 




(P/Po) 


(T/T^f (M/Mj^(p/pJ^ 


120 



Figure 45, and the analysis of Reference 8, that the boundary layer displace- 
ment thickness can be determined by a real gas analysis and thus the problem 
is tractable. In order to get numerical results, the Re 3 molds number, Rey^, 

and the Mach number, shoAvn in Figure 45 must be recogmzed as those 

correspondii^ to the unsubscripted parameters M, p, T, etc. in the analysis 
of Table Vni (i. e. , the local values downstream of the obliq,ue shock) . There- 
fore, the solution will, of necessity, be an iterative one similar to that 
described in the boundary layer -shock interaction analysis of Reference 2. It 
IS clear from Figures 56 through 59 that this refinement in the shock wave 
interaction problem is essential if accurate aerodynamic heating rates are to 
be calculated. The correlation of pressure rise and aerodynamic heating can 
still be obtained in the form shown in Figures 51 and 52, 


Conclusion 


From this discussion in which the various factors involved in the 
determination of the aerodynamic heating of bodies traveling at hypersonic 
speeds are presented, it should be evident that even the "simple" mathe- 
matical analyses and physical models rec[uired to get an approximate solution 
are quite complex and can only be found by numerical techniques Although 
more accurate analyses and physical models can be used at each step in the 
calculation (inviscid real gas solution, the boundary layer analysis, and 
the real gas interaction problem, etc.) the questions of computer capacity 
and cost versus the desired engineering accuracy must be considered. 
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FIGURE 42. COMPOSITE PLOT OF VELOCITY PROFILES WITH T /T = 1/3 

FOR AN IONIZED MONATOMIC GAS ^ 
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FIGURE 43. COMPARISON OF FROZEN AND EQUILIBRIUM TEMPERATURE PROFILES 
WITH T^/T^= i/3, 0, e^= 0.997, and M^= 5 FOR AN IONIZED MONATOMIC GAS 
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FIGURE 44. COMPOSITE PLOT OF TEMPERATURE PROFILES WITH T /t^= i/S FOR AN 

IONIZED MONATOMIC GAS ( Frozen Flow) ^ 
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FIGURE 46. VARIATION OF e WITH RESPECT TO u/V WITH T /T = l/S, 
€ = 0, € = 0.600, AND M = 5 FOR AN IONIZED MONATOMIC GAS 
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FIGURE 47 IONIZATION PROFILES WITH T /t = l/s, e = 0, e = 0 600 AND M 

W °° W *0 

FOR AN IONIZED MONATOMIC GAS 
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FIGURE 48. EFFECT OF BOUNDARY-LAYER THICKNESS ON THE DISTRIBUTION OF h/h WITHIN 2 l/2 

o 

INCHES OF THE FLAT-PLATE CENTER LENE FOR A 2 . 8-INCH-DlAMETER RIGHT CIRCULAR CYLINDER; 

M = 3.51; R « 2.90 x lO®. 


130 


BOW 



FLAT PLATE 


FIGURE 49. MAXIMUM HEATING — BOUNDARY LAYER-SHOCK WAVE INTERACTION 
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5IGURE 50. SHOCK INTERACTION WITH TURBULENT BOUNDARY LAYER 
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FIGUEE 51. PRESSUEE RISE VERSUS MACH NUMBER - REYNOLDS NUMBER 
PRODUCT FOR LAMBDA-FOOTED SHOCK 
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o experiment Ref. I 
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FIGURE 52. MAXIMUM HEATING RATE AS A FUNCTION OF 




pressure rise ratio, p/p 
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FIGURE 54. PRESSURE DISTRIBUTION IN LAMINAR AND TURBULENT SEPARATED FLOW 
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FIGURE 57. OBLIQUE SHOCK TEMPERATURE RATIO INCLUDING 

REAL GAS EFFECTS 
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Introduction 


Thermophysics research at Marshall Space Flight Center is concerned 
with the thermal space environment and its effect on space vehicles. Space 
craft are in radiative equilibrium with their environment, the main factors 
being solar radiation, the planetary albedo, the planetary infrared radiation, 
and aerodynamic heating only during the final phases of a decaying orbit. The 
following paragraphs deal with various aspects of thermophysics research. 

The emissivities of various surfaces of space vehicles are essential 
to control the radiative exchange with the environment of celestial bodies like 
sun, planets, or moon. 

Theoretical studies enable us to est^lish mathematical models and to 
develop computer programs for thermal design and for environmental effects, 
to other aspects of thermophysics. 

Research in radiative and conductive physics is an experimental 
program carried out to solve such problems as modes of heat transfer m 
underdense materials, interfacial conductance between materials in space and 
others. 


Thermal control has become a vital design area for all space vehicles. 
My division has been fully responsible for the thermal control of Explorers I 
through V, VII through XT, upper stages of the Saturn I flights, SA-4 and SA-5, 
and for the Pegasus satellites I, n, III, All of our research finds application 
in the thermal control of such spacecraft and contributes in general to the 
advancement of technology in this vital area of space exploration. The area 
includes such problems as the determination of launch windows, the selection 
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of thermal control coatings inside and on the outer skin of spacecraft and 
computer programs to predict temperatures for the total lifetime of the 
satellite in space. 

The effects of the space environment on thermal coatings proved to be 
a very interesting research area with many problems still unsolved today. 

Even if we have found a good thermal design, if we have verified it in thermal- 
vacuum tests, and have predicted the changes of temperature in space due to 
orbital precessions, earth motion with respect to the sun, etc. , we have to 
ascertain that the initial values of material emissivities do not change in space. 
We found that this is not possible Even the most stable coatings and paints 
that we developed undergo considerable changes. The research emphasized 
two aspects: (1) to understand the basic degradation mechanisms, and (2) to 
malce the changes due to the space environment predictable, such that it can be 
taken into consideration by the thermal design engineer. 

An important aspect of direct involvement in the thermal control of 
spacecraft is the feedback of knowledge from the actual performance of the 
spacecraft in orbit. The evaluation of temperature measurements helps to 
learn how good the design approach was, how close the predicted temperature 
fluctuations are matched by the actual values. However, it proved to be very 
difficult to pinpoint the reasons for deviations, because more than 20 variable 
parameters are involved. We, therefore, initiated the development of a flight 
experiment of "thermal environment sensors. " Various versions of the experi- 
ment were flown on the early Explorers, on the Saturn upper stages and on all 
three Pegasus satellites. Results from these experiments contributed to 
explain discrepancies between lab experiments and actual exposure to the space 
environment 

Because of the great demands for better knowledge about emittances and 
reflectances of solid matter for space application, new and advanced methods 
for laboratory measurements had to be developed. These include IR 
spectroradiometers, hohlraum reflectometers, calorimeters, in situ 
irradiation with UV and simulated solar wind and simultaneous measurement 
of the optical characteristics, bidirectional devices which allow to have 
goniometric variations of the incident and reflected beams. 

Table IX gives a summary of the "work areas" in thermophysies 
research which are pursued by the Space Thermophysics Division. More 
details will be given in subsequent chapters . Activities in the "work areas" 
cover a wide range from basic research to the application for MSEC mainline 
space flight projects. The area of space flight experiments has been an 
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especially rewarding effort, because of the stimulus on our laboratory research, 
our theoretical studies, and its direct tie-in with the scientific payloads The 
continuous participation in the Marshall mainstream efforts and the application 
of our research to space flight projects has proven to be rewardii^ in a dual 
sense- (1) the contribution to the MSFC projects makes our work an integral 
part of these projects, and provides an excellent training for our physicists 
and engineers, and (2) the experience gained allows us to plan future research 
and the selection of scientific areas with potential for contributions to the space 
flight effort. 
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TABLE IX. THERMOPHYSICS DIVISION WORK AREAS (Space Sciences Laboratory, MSEC) 


PROGRAMS 

ACTIVITIES 

Thermal Similitude 

Mathematical studies of new approaches to thermal 
model testing 

Theoretical Studies of Lunar Thermophysics 

Theory of electromagnetic waves interacting with the 
lunar surface, phase relationships and amplitudes of 
lunar subsurface heat waves 
Evaluation of IR and thermal mm wave measurements 
of the moon 

Electromagnetic Radiation Physics 

Especially IR and far IR 
IR of the moon 

Thermophysical Properties of Underdense Materials 

Thermal, electrical and acoustic contact resistance 
between particles 
Reduced gravity environments 
Conductive and radiative physics 

Lunar Geological and Geophysical Studies 

Mechanical forces between particles of simulated lunar 
materials 
Lunar drill studies 

Emissivity Physics 

Advanced measuring techniques of surface optical 
properties 

Effects of micrometeoroids 
UV and solar wind effects 




table IX. THERMOPHYSICS DIVISION WORK AREAS ( Concluded) 


PROGRAMS ACTIVITIES 

Solid State Physics and Reaction Kinetics Interaction of photons with solid matter 

Degradation mechanisms of thermal control coatings 
EPR and NPR measurements of unpaired electrons 
Magnetic susceptibility studies 
Solar UV measurements 

Thermal Space Environment Solar simulation 

Bidirectional albedo studies 

Theory of Rayleigh scattering and its effect on space 
vehicle thermal control 
Thermal control of space vehicles 









Thermal Similitude 


The thermal vacuum testing of space hardware requires larger and 
larger chambers as the sizes of spacecraft increase. It is, therefore, of 
considerable interest to establish nondimensional modelii^ laws which allow to 
do some of the thermal testing with models of space hardware. During the 
past several years, research has at many places in the country unfortunately 
focussed on the question whether the models can replace the full scale tests. 
Some of the results seemed to disprove the usefulness of thermal models, 
however, their results are not conclusive. Research at the Space Sciences 
Laboratory was aimed at establishing dimensional parameters for model 
testing in a similar way as is done in wind tunnel tests for aerodynamic 
profiles. Just as the wind tunnel test does not replace the full scale test flight 
but helps the designer with vital knowledge about aerodynamic coefficients, 
thermal model testing can be expected to help in the early phase of the space- 
craft design to establish the thermal design concept. A computer program was 
set up to find non-redundant and consistent sets of tt - factors which are the 
nondimensional quantities . 


For thermal modeling four tt ratios and for time variable systems five 
7T ratios have to be used. The computer program yielded a total of 57 different 
sets of such IT ratios. The following equations show such a a set of tt ratios 
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The symbols are defined as follows: 


T = temperature of Jth node, [©] 

3 

= heat capacity of jth node, [H0 



over-all conduction coefficient between nodes k and j 
[HT“^ ©"^] 
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J 
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t = time, [T] 

= area of jth node radiating to space, [L^] 
j t 

■- temperature of kth node, [ O] 

The symbols for the dimensions are defined as follows: 


overall radiation coefficient for net radiative heat transfer 
from node k to j, [HT ^ 6 

“1 -1 

thermal input intensity on the jth node, [HT L ] 
area of jth node receiving thermal input, [L^] 


© = temperature in 

H = heat or energy in joule 

T = time in seconds 

L = length in meters 


We assume 


7T as a function of the other tt ratios, 
s 


That is, 


TTg - f ( TTi, 1T2 j 7T3, 1T4 ) . 


If we solve for T. then. 
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We have concentrated on transient thermal problems because of their impor- 
tance for near earth orbits. The ingress into and egress out of the earth’s 
shadow calls for model investigations of thermal transients. The theoretical 
studies were supplemented by experimental work on a set of models with well 
defined but simple geometries Besides the largest model or "prototype" 
two scaled down models with linear sealing factors of 1 : 2 and 1 • 4 were used. 
Figure 60 shows the arrangement of a conical model. The top and bottom 
plates of the cone are connected by four rods which support an additional mass 
between them. The figure also explains the numbering system of the thermo- 
couples Figure 6 i shows the cyclic response of the small top plate of the 
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FIGURE 60, CONICAL MODEL USED IN THERMAL SIMILTUDE TESTS 


coti© frustum (thermocouple #i) due to the imposed heat indicated by the 
arrows in Figure 61. Figure 62 shows the thermal transients of the block 
supported between the four rods.' A built-in heater is able to generate very 
steep temperature rises. The more graduate drop of temperatures is due to 
radiative and conductive heat transfer. The temperature curves of the proto- 
type and the two scaled-down models follow the same pattern The absolute 
match of the temperatures is not essential. However, the modeling laws should 
allow to predict the temperature transient for each, if the initial boundary 
conditions are given. The research has revealed that thermal models can be 
helpful to solve specific questions. Since the testing for full-size space vehicles 
in vacuum chanbers is quite extended in time and costly, not all necessary 
parameters and thermal transients can be tested However, they can be 
tested in small models for a fraction of the cost and time, and as such 
supplement the full-scale tests. 


Theoretical Studies Of Lunar Thermophysics 


Theoretical studies are carried out to establish a better understandmg 
of the interaction of electromagnetic waves with the lunar surface. Such an 
understanding is necessary to determine the environmental criteria for lunar 
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FIGURE 6i CONFIGURATION C ^ TRANSIENT THERMAL RESPONSE OF 

THERMOCOUPLE NO. 1 

surface craft. The lunar surface has very unique properties which are quite 
different from those found for natural surfaces on the earth. One of the 
properties is the strong directivity of the backscattered solar radiation or 
photometric function. Figure 63 shows this function as derived by Hapke. It 
very closely approaches the reflected solar radiation for an average element 
on the lunar surface. It is a photograph of a.^three-dimensional model. The 
downward arrow is the direction of the sun's incoming radiation at an angle 
of 30° from the local vertical. The upward pointed arrow gives the local 
vertical. The intensity of the scattered or reflected radiation is depicted by 
the length of each point of the three-dimensional figure from the intersection of 
the sun vector with the local vertical in the center of the horizontal plane. 
Hapke's function is given by the equation also shown on Figure 63. 
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TIME-MINUTES 

FIGURE 62. CONFIGURATION C - TRANSIENT THERMAL RESPONSE OF 

THERMOCOUPLE NO. 3 
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FIGURE 63. GEOMETRICAL REPRESENTATION OF THE HAPKE FUNCTION 
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B(a, g) = 1 for a 
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The angles in equation (1) are explained in Figure 64. The same 
figure shows an equation for the infrared emission flux (radiance) from a lunar 
surface element as a funccion of the sun elevation ar^le i, observer elevation 
angle e and the lunar phase angle a. The angle is defined by the second 
equation as function of the known angles i, e (p^, <p^. A^, A 2 , A 3 are empirically 

adjusted constants. The third equation of Figure 64 is the Stephan-Boltzmann 
equation which allows to compute radiance temperatures from the radiative 
Flux B, as can be seen from Figure 65. This radiance temperature is 
dependent on the angle of observation. It shows this temperature for an 
incidence angle i of 30® for an azimuth angle of the observer of zero, i.e. , 
co-planar with the sun vector. The concept of a temperature which varies with 
the observation mi^t be very disturbing. If we could assume that the tempera- 
ture of the lunar surface spot is well defined and uniform, we could introduce 
such a temperature and compute the emittance which is variable with the 
viewing angle. However, this is not possible at the present time, because we 
still have some serious open questions which do not allow to define a constant 
temperature. One of the questions is whether the radiative behavior is due to 
the lunar structure on a microscale as inferred in the Ashley expression shown 
on Figure 64 or whether it is due to shadow effects in craters on a macroscopic 
scale. Another question is whether a meaningful lunar surface temperature can 
be defined, if we have a very underdense material on the surface with a 
temperature gradient of the order of 100®K/cm. Figure 66 shows the brightness 
or radiance temperatures as a function of the elevation angle of the observer 
with various sun angles as parameters. It can be seen that the brightness 
temperature peaks in all cases close to the direction of the sun angle. The 
theoretical curves are compared to measurements taken by various authors 
which are referenced on the graph. 

Figure 67 shows a plot of the temperature of the full moon across the 
disk. Various sets of measurements are shown together with the Lambert 

distribution of temperature which should follow a cosV '*0 relationship and the 

empirical relationship of a cos^'^®^ law established by Pettit and Nicholson. 

The deviation from Lambert's law has been explained that it is due to the 
surface roughness of the moon. 
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FIGURE 64 ANGULAR EXPRESSIONS USED IN THE HAPKE FUNCTION 

Th.e problem of the variation of lunar surface temperatures during a 
lunation is important for the planmng of lunar explorations. A computer 
program has been developed to compute temperatures using the thermal 

parameter (kpc)”^/^. Figure 68 shows the result of such computations for 
the landing site of surveyor I. Also included in the graph are curves based 
on measurements by various, investigators . The measurements from Surveyor 
show strong deviations. This has not been explained yet. It is quite possible 
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that the stror^ angular dependence and the above mentioned prohlems in connec- 
tion with the definition of lunar temperature have a bearing on this discrepancy 
between measurements and our understanding of the moon expressed m the 
computed curves. It is at this moment not possible to match the computed 
curves by simply varying the values of the thermo physical parameters. 



MGURE 65. LUNAR BRIGHTNESS TEMPERATURE AS FUNCTION OF THE 
ANGLE OF OBSERVATION FOR SOLAR INCIDENCE OF 30“ 

Electromagnetic Radiation Physics 


The theoretical studies described previously are closely related to the 
experimental results obtained from ground based measurements made by 
various observatories. We recognized early that simultaneously with a 
theoretical effort we had to simulate the conditions in the laboratory, and to 
make additional observations of lunar IR radiation with ground based telescopes. 
Our Electromagnetic Radiation Laboratory covers a broad spectrum from solar 
radiation to the thermal infrared (IR) through the far IR extending all the way 
to the mm wave region. The instrumentation includes various IR methods such 
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as IR radiometry, grating, and filter spectrometers, Miclielson interfero- 
metery, and IR imagery. Another area of instrumentation is that of IR 
detectors . Our investigations cover a wide variety of quantum and bolometer 
detectors which are needed for the various spectral bands of IR. An example 
is the gallium doped germamum bolometer detector which is suitkble for the 
far IR from wave members of 200 to 5 cm ^ I would like to show some 
examples of the instrumentation. Figure 69 shows the schematic of an IR 
radiometer which uses an IR filter with a band width from 10.2 to 12jLtm. 
Figure 70 shows this radiometer mounted on a Cassegrainian telescope The 
research program is carried out in connection with extraterrestrial resources 
to determine the IR and thermal properties of lunar and planetary soil and 
minerals . Figure 71 shows an IR scan over the crater Tycho during the 
penumbral phase of the lunar eclipse of April 12, 1968 Such rayed craters 
as Tycho are showing as hot areas against a cooler surrounding shortly after 



FIGURE 66. BRIGHTNESS TEMPERATURE VERSUS OBSERVER 

ELEVATION ANGLE 
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sunset during a normal lunation cycle, but even more pronounced during an 
eclipse Figure 72 shows a schematic of a Michelson interferometer covering 

the spectral range up to 2000 jLun. The 
advantage of this instrument is in its 
integration times shorter than any other 
instrument in order to obtain a complete 
spectrum. In grating instruments most 
of the energy is eliminated except for 
one spectral line or band, whereas the 
interferometer allows to handle the 
complete spectrum simultaneously 
Figure 73 shows a UBV photometer 
which is used for color photometry of 
stars. The letters TJ, B, and V stand 
respectively for the ultraviolet, blue, 
and visible filters, whose spectral 
characteristics are shown in Figure 74 
together with those of other filters 
extending into the IR The scientific 
project is concerned with the measure- 
ment of binary stars and the analysis of their orbital elements The UBV data 
will be supplemented by a photometer which operates in the various IE. bands 
shown in Figure 74. This is specifically important for cases where the 
comparison star is a late type star with a low temperature which corresponds 
to a blackbody maximum in the IR. Figure 75 shows some results obtained at 
the University of Pennsylvania with IR photometry It can be seen that the 
secondary minimum of the eclipsing binary almost disappears m the noise for 
visible light, but is quite pronounced in the IR. 

Figure 65 discussed previously has been obtained by ground-based 
measurements. Similar results were obtained by goniometric IR measurements 
in the Electromagnetic Physic s Lab Figure 76 shows the experimental 
instrument in the laboratory used for these measurements 

Thermophysical Properties Of Underdense Materials 
And Interface Conductance 



FIGURE 67. DISTRIBUTION OF 
PLANETARY HEAT ACROSS THE DISK 
OF THE FULL MOON 


Two problems connected with space vehicle projects at MSFC made the 
investigation of underdense materials attractive: Insulating powders which 
could accomplish required reduction of heat transfer under difficult geometrical 
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FIGURE 68 . CALCULATED AND MEASURED TEMPERATURES FOR 

SURVEYOR I SITE 


situations and the investigation of the top layer of the lunar surface. The 
properties of insulating powders have to be known under zero-g conditions, 
whereas lunar conditions correspond to i/6-g. Our studies started initially 
with the problems of insulation materials for thermal design application. 
However, with the initiation of the lunar exploration by Ranger, Surveyor, and 
expecially project Apollo the thermo physical properties of the moon took hi^er 
priority. Another aspect was that our thermal design studies showed that the 
insulation properties could be achieved by 'well known means other than 
powdered insulation materials . In connection with the lunar exploration it 
seemed desirable to simulate expected lunar surface materials in the laboratory 
and study the possible effects of reduced gravity by analytical means. The 
thermo physical properties of the moon are needed for correlation with the 
mechanical properties of the soil and also for correlation with the radiative 
characteristics of the moon m the IR and in thermal mm waves, which already 
has been discussed. Our research concentrates on the study of the heat transfer 
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RAD I ATI OH FROM TELESCOPE 



FIGURE 69 . INFRARED RADIOMETER 

modes, on their temperature dependence, on the influence of vacuum grain 
size, type of material and density. Figure 77 shows the dependence of the heat 
transfer on the density. At high densities the conductive transfer is dominant, 
it decreases with decreasing density. For hi^ly underdense materials the 
conductive term is nearly negligible, but the radiative transfer is hi^. In 
between is a minimum Present thinking about the lunar surface is that a thin 
layer on the top is highly underdense and that the packing density increases 
with depth. This is due to the ''gardening effect" of the continuous impact of 
meteoroids and micrometeoroids After one or several meters the density 
approaches that of solid rock Research in this area is carried out in our 
Conductive and Radiative Physics Lab. 

Another research activity conducted in this laboratoiy is in connect- 
ion with the interface resistance between metallic surfaces under the hard 
vacuum conditions of space. The so-called metallic contact between two planar, 
well polished surfaces is in reality due to the airfilm between both surfaces. 

If this film is removed, the heat transfer depends on the conduction through a 
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HGURE 70 . LUNAR IR EMSSION DATA ACQUISITION USING AN IR 
RADIOMETER ON UNIVERSITY OF GEORGIA 24” CASSEGRANIAN TELESCOPE 

few actual grain to grain contacts and on radiative transfer \\*ich is low for 
polished metals. The interface resistance could become a real problem for 
metal to metal contacts which have to carry a heat load in space. Figure 78 
shows some results of heat transfer of an interface as a function of the pressure 
applied. The research is presently concentrated on a better theoretical under- 
standing of the surface physics problems involved. 


Lunar Geological And Geophysical Studies 


Project Apollo has the purpose to land two astronauts on the moon and 
to bring them back safely. The problems connected with the landing require 
a knowledge of the mechanical properties of the lunar surface and the 
geological features and composition of the moon. A laboratory for geological 
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FIGURE 71. INFRARED SCAN ACROSS ENTIRE LUNAR DISK DURING 

UMBRAL PHASE OF ECLIPSE 

and geophysical studies was set up for such studies. Figure 79 shows an 
ultraclean vacuum chamber which can be evacuated to 10 ~ torr. It provides 
for effective outgassing of simulated lunar ^^rface material by tumbling of a 
drum. The vacuum was maintained at 10 torr even during such an out- 
gassing experiment. The materials were recommended by Dr. E. Schoemaker 
of the U. S. Geological Survey and other geologists who study lunar properties 
such as finely ground powders of basalt, pumice, ryolite, etc. The results 
were very interesting . Basalt powder starts to adhere to the steel wall of the 
chamber and the grains to each other shortly after the outgassing is started. 
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HGURE 72 MECHELSON INTERFEROMETER 





FIGURE 73 COMPLETELY DIGITIZED MULTI-COLOR PHOTOMETER 
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FIGURE 74. FILTER SPECTRAL BAND PASS AND LETTER DESIGNATION 

OF JOHNSON SYSTEM 


After about 15 minutes, all of the powder adheres to the wall and does not drop 
anymore This explains the consistency of lunar soil to behave somewhat like 
wet sand or the soil of a freshly plowed field. Since there is no amount of free 
water on the moon, the adhesion could not be explained in terms of capillary 
forces and adhesion of water to the soil grains. However, if the lubricating 
and separating surface layer of air is removed from silicate type minerals or 
from any other material, the van der Waals forces, chemical binding forces or 
electronic forces of metallic bonds become effective. It is this effect which 
makes the lunar "dust” behave more like "soil.” Results found in our 
laboratories were continued by and are in agreement with the results of the 
recent Surveyor experiment of the Soil Mechamcs Surface Sampler. 

Another example of our research in this area is the analysis of 
simulated lunar materials by X-ray diffraction and Bragg crystal techniques. 
This area is quite promising to provide support for the X-ray solar astronomy 
projects for the Apollo Telescope Mount. Members of this Laboratory are 
the experiment scientist and alternate experiment scientist on S-054 which has 
Dr. R. Giacconi of the American Science and Ei^ineering Corporation as 
Principle Investigator 
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Emissivity Physics 


The selection of the emittances of outside surfaces of the spacecraft is 
a powerful method to control the' onboard temperatures. The area of thermo- 
phyiscs obtained a considerable impetus because of this challengii^ task. In 
addition to the development of space-stable thermal control surfaces for a 
variety of different requirements, members of my division were involved in 
research in emissivity physics. New and advanced methods of laboratory 
apparatus were developed and our experimental facilities include a broad 
capability including hohlraum apparatus, radiative calorimeters, radiometers, 
and reflectometers covering ultra violet, visible and infrared radiation. This 
paper describes one of the devices, a bidirectional spectroreflectometer 
(Fig. 80) . The arrows indicate the rotational degrees of freedom. The 
sample holder is mounted horizontal and allows investigation of any type of 
materials including powders such as underdense material simulating the top 
layer of the lunar soil. The sample can be rotated around its vertical axis and 
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is illuminated with monochromatic li^t obtained from a prism monochrometer 
or from a laser. The illumination angle can be varied in one plane with nearly 
a complete hemicircle (±90 degrees) . The measuring head can be placed 

anywhere within a hemisphere by 
adjustment of the azimuth angle from 0 
to 360 degrees and the elevations from 
0 to 90 degrees. The backscattering 
direction which is usually not included 
in such equipment is, in our case, 
obtained by a beamsplitter arrange- 
ment. Precautions have to be taken that 
the beam splitter does not pick up 
radiation from the surrounding at 90 
, degrees to the incoming illuminating 
radiation. Figures 81, 82, 83 give 
results for MgO with the incident beam 
and the measuring beam in one place . 

All measurements are for infrared 
radiation at \ = 2. 0 pm. The three 
figures differ in the incidence angle, 
which is 0 degree in Figure 81, -30 
degrees in Figure 82 and -60 degrees 
in Figure 83. Magnesium oxide is an 
especially interesting material because 
it is used as a standard and as diffusing 
surface material for integrating spheres. 
It can be seen that the reflection of MgO 
follows closely that of a Lambertian 
surface for all reflection angles with the 
exception of the backscattering angle. 
Because it is difficult to measure at 
the backscattering angle, this effect had not been discovered. It is interesting 
to note that the curve ends in a sharp peak at the exact backscattering angle. 

This also resolves a problem in connection with the photometric function of 
the moon. Because the earth is in the way, the closest angles that can be 
measured are between 1.5 to 5 degrees before and after lunar eclipses. It was 
assumed by many investigators that the reflected li^t follows a shallow 
maximum and therefore the curves were rounded off at the backscattering point, 
instead of ending in a sharp point. This question has been resolved in favor of 
a sharp peak at a 0 -degree angle between incident and measuring beam. 



FIGURE 79. VACUUM CHAMBER 
USED IN STUDIES OF POSSEB LE 
LUNAR MATERIALS 


170 



FIGURE 80. BIDIRECTIONAL SPECTROREFLECTOMETER 


171 


iA 

Z 

D 


a: 

< 

c 

K 

s 


D 

0. 

K 

3 

O 

q: 

o 

K 

U 

UJ 

H 

UJ 

o 



-60 -50 -40 -30 -20 -10 0 10 20 30 40 50 60 

B (DEGREES) 


FIGURE 81. REFLECTANCE OF MgO VERSUS ANGLE FOR ILLUMINATION 

ANGLE = 0 

Another example of the activity in emissivity physics may be given; 
the measuring of optical properties of thermal control coatings prior to and 
after exposure to ultraviolet irradiation. Figure 84 shows the spectral 
reflectance of a zinc oxide-silicone paint for the wavelength range from 0, 3 to 
2.4 iim. Measuring points and curves are shown for initial conditions and for 
various time integrals of equivalent sun hours obtained by the irradiation of 
the sample with an A-H6 UV lamp. All curves show the absorption edge in 
the UV end of the spectrum due to the band gap of 3.2 eV of ZnO. The 
degradation by the UV irradiation is most pronounced in the red and IR end of 
the spectrum . The measurements were made with a spectroreflectometer . 
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FIGURE 82. REFLECTANCE OF MgO VERSUS ANGLE FOR ILLUMINATION 

ANGLE = 30° 


Solid State Physics And Ultraviolet Physics 


Our concern with the space environmental effects on thermal control 
coatings necessitated an investigation of the basic mechanisms that underlie 
such effects as the degradation (yellowing) of "white” thermal control surfaces 
During the development of the early Explorer satellite I had already a healthy 
respect for the environment. Very little was known at that time and surfaces, 
such as sandblasted metal and AI 2 O 3 , which were assumed to be very stable, 
were selected for the outside skin of Explorers I through XI. When satellites 
became larger and more sophisticated, new coatings were needed. Some of 
the most successful coatings developed by one of our contractors, the Research 
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FIGURE 83. REFLECTANCE OF MgO VERSUS ANGLE FOR ILLUMINATION 
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Institute of IIT, were S-13 (which has ZnO as pigment and a methylsilicone 

as binder) , S-13 G (an improved version of S-13 and a silicate film on the 

ZnO grains) , and Z-93 (which consists of SnO and silicates) . Some essential 

requirements were easy applicability to large surfaces and to be as space stable 

as possible with respect to UV radiation, temperature, and vacuum. The 

stability with respect to solar wind becomes important for deep space probes 

or very eccentric satellites. In many cases coatings with best stability for 

UV radiation are still affected by solar wind and vice-versa. Actually, none of 

the coatings showed completely satisfactory behavior. The white S-13 paint 

on the S-IVB stage of the Pegasus satellite had an initial oi /e. value of 0.22 

s ir 

which increased during the launch trajectory and satellite injection into orbit 
to 0.51. This large increase caused the stage to have a much higher 
temperature than intended for the thermal control of the satellite.^ 

We have been instrumental in the formulation of many of the nation- 
wide research objectives and have contributed with original research in our 
.laboratory. Some of the laboratory instruments used are showm in the 
following figures. Figures 85 and 86 show the 2 electron parametric resonance 
(EPR) apparatus operating in the X-band (8 to 12 x 10^ Hz) and K-band 
(3. 5 X 10^° Hz) . Both have a double cavity setup, one for comparison and one 
\idiich can be irradiated with UV v^ile the measurements are performed. Figure 
87 shows some results obtained with the X-Band EPR. The frequencies are 
swept by varying the magnetic field H. The resonance is due to UV irradiation 
which generates unpaired electrons. Other instrumentation includes an 
apparatus for magnetic susceptibility. Hall effect equipment, UV irradiation 
apparatus for thermal control samples, solar wind simulator, thin film 
measurements, and photoconductance apparatus. 

Another aspect of our activity is the definition of the UV radiation 
environment for space vehicles and the participation as experiment scientists 
in the implementation of present or future scientific experiments concerned 
with the measurements of solar UV and stellar UV. Figure 88 shows a UV 
spectrum of the sun compared to a spectrum measured by the Naval Research 
Laboratory (NRL) with the Zeta facility (Harwell) . Figure 89 shows a 
spectrograph by the NRL which covers the spectral region from 500 to 1500 A. 

A schematic of the S-083 B solar astronomy experiment of NRL on the Apollo 
Telescope Mount A(ATM-A) , is shown in Figure 90. MSFC is responsible for 


1. An overheating of the satellite and especially of the electronics canister was 
avoided because of an active louver control system which was introduced 
into the thermal control concept as insurance against such unforeseen 
occurrences. 
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FIGURE 85. CLOSE-UP OF MAGNETS OF ELECTRON SPIN 
RESONANCE EQUIPMENT 



FIGURE 86. ELECTRON SPIN RESONANCE DATA-ACQUISITION 

EQUIPMENT OF SSL 

the implementation of the ATM, the first manned astronomy project of the 
Apollo Applications Program. My division is responsible to furnish experiment 
scientists who work with the Principle Investigators ( P. I. ) of the NRL, 

Dewitt Purcell, and Richard Tousey. An experiment which is being prepared by 
another P.I. of the NRL is in the early stages of development. MSFC is 
presently studying this and many other experiments approved by the Astronomy 
Subcommittee of the Office of Space Sciences and Applications (OSSA) . Figure 91 
shows a schematic of the experiment. It is a double reflecting instrument 
planned for an astronomical sl^ survey in the far UV region from 900 to 
1800A. The low efficiency of the reflection grating and the primary mirror are 
compensated by the use of an image converter. The secondary mirror of the 
Cassegrain telescope is the photocathode. The electrons are accelerated by 
20 kV and form an image on the film in the Cassegrainian focus with the help 
of a focussing magnet. Figure 92 shows a photograph of the assembled 
instrument. It was flown suecessfully in a rocket flight at White Sands Proving 
Ground, and has proven the principle of operation and furnished interesting 
far UV results of several stellar sources. 
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Thermal Space Environment 


This research area covers one of the main assignments of the Space 
Thermophysics Division. The definition of the space environment is of 
importance to all spacecraft and affects their design and operation. The 
simulation of the thermal environment in the laboratory and the study of thermal 
control will also be described in this connection. 
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THE SPECTRUM OF THE SUN (US NAVAL RESEARCH LABORATORY) 
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FIGURE 88. XUV SPECTRUM OF THE SUN FROM 160 TO 220 A 
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ASTIGMATISM) 

THE DOUBLE -DISPERSION SPECTROGRAPH OF DETWILER, 
GARRET, PURCELL, AND TOUSEY (1961), FIRST FLOWN ON 
APRIL 19,1960. 

FIGURE 89. DOUBLE DISPERSION SPECTROGRAPH OF NRL 


179 








FIGURE 90. XUV SPECTROGRAPH, NRL ATM EXPERIMENT B 

Figure 93 shows the result of computations of Rayleigh scattering in the 
earth's atmosphere, and the irradiation of a satellite surface element by the 
earth albedo due to Rayleigh scattering. The figure shows that the spectrum 
of the scattered sunlight differs from the solar spectrum. It is interesting to 
note that the spectral albedo radiation depends on the inclination of the 
satellite surface to the sun vector and to the radius vector to the earth. Figure 
94 shows a radiometer for solar, albedo, and earth IR radiation. This concept 
has been studied as possible contender for a flight instrument. The solar 
constant has, to date, only been determined by extrapolation of ground 


180 




FIGURE 91. SCHMIDT IMAGE CONVERTER STELLAR SPECTROGRAPH 

OF NRL 


measurements to space. The specific device shown here will not be selected; 
however, it shows some of the basic features of such radiometers: a black 
body cavity, insulation from the housing, thermal guards, and the possibility 
to heat with measured amounts of electric energy. One of the requirements 
for solar radiometers is a controlled attitude. Its control does not have to be 
very ti^q however, the attitude has to be known to an accuracy of at least 
a degree. This allows. together with the ephemeris of the satellite orbit, to 
'compute separately the solar constant, albedq and IR radiation. 

Figure 95 shows a solar simulator. Its cylindrical vacuum chamber 
can be evacuated to 10 torr and the walls can be cooled to liquid nitrogen 
temperature. An Aerospace Controls Corporation xenom-type solar lamp allows 
simulation between 0.5 and 2 solar constants. The solar simulator is used for 
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FIGURE 92. FUGHT PACKAGE OF SCHMIDT SPECTROGRAPH (NRL) 
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FIGURE 93. RAYLEIGH SCATTERING IN THE EARTH'S ATMOSPHERE 
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DETECTOR ASSEMBLY THE BLACKBODY IS IMBEDDED IN A THERMOSTAT SET 
TO THE MIDRANGE OF EXPECTED TEMPERATURES 

FIGURE 94. RADIOMETER USED FOR SOLAR ALBEDO AND 

EARTH RADIATION 

investigations of the radiative equilibrium of space components and the testing 
of thermal shields such as the frontal area of a solar telescope for simulation 
of the irradiation of simulated lunar materials through a lunation or through 
a lunar eclipse . 

Figure 96 shows one of the three Pegasus satellites. All three were 
launched in 1966. My division had the responsibility for the thermal control of 
the Pegasus experiments and of the electronics canister. In addition to the 
primary engineering experiments about micrometeoroid penetration 
frequencies, the three satellites carried thermal experiments with four sensors 
each for determination of environmental effects on thermal control surfaces 
(Fig. 97) . The location of the thermal environmental effects experiment is 
indicated on Figure 96, The sample surfaces varied from satellite to 
satellite and included, in all cases, two of the coatings which were used for 
thermal control of the Pegasus satellites. These were the S-13 and the Alodine 
thermal control surfaces. S-13, which was used on the Saturn S-IVB stage, 
is a white coating which is highly reflective for solar radiation and contains a 
zinc oxide coating with a silicone binder. The Alodine coating is a greenish 
electro-deposited coating used on the wings of the Pegasus. 


184 




FIGURE 95. SOLAR SIMULATOR AT SSL 

A new concept of preflight measurements of satellite thermal control 
coatings was introduced with the Pegasus satellites. The time interval between 
the application and vacuum testing of thermal control coatings and the actual 
launching is usually 4 to 6 months but may extend to more than a year. No 
provisions are made to reconfirm or test the accuracy of the emittance values 
of the satellite thermal control surfaces. For the Pegasus satellites we 
introduced a portable device which allows measurement, under "-field 
conditions," of the solar reflectance or the IR emittance with two different 
measuring heads. Figure 98 shows the solar reflectometer head being used on 
flight hardware. Measurements with the portable instruments were made on the 
Pegasus satellite at various sites such as MSFC, Fairchild, the satellite prime 
contractor, and finally in the hangar and on the launching pad up to two days 
prior to the launch at the Kennedy Space Center. The status of all thermal 
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FIGURE 96. PEGASUS SATELLITE 



FIGURE 97. THERMAL EN\HRONMENT SENSORS FLOW^ ON PEGASUS 
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control coatings prior to satellite launching could thus be ascertained, and this 
information was used for comparison to on-board measurements of the 
degradation by the space environment. We know from this that the initial rapid 
decrease of solar reflectance is a true space effect or atmospheric environ- 
ment effect during launch. 

A promising approach for the thermal control of easily overheated 
components of experiments on a boom-like magnetometer is the use of fusible 
materials. Requirements are a high heat of fusion and a minimum amount of 
supercooling effects for the fusing material. Figure 99 shows a flight experi- 
ment which was assigned the NASA Headquarters nomenclature Number T005. 

It is a schematic of the adiabatic test model of the heat of fusion radiator. The 
fusion material is embedded between the meshes of a honecomb system for 
increased heat conduction. Figure 100 shows several temperature -time curves 
for different rates of heating per unit area. The main advantage of the heat-of- 
fiision principle is that it can be used in addition to other passive or active 
thermal control concepts, and that it reduces the extremes in temperatures. 
For this reason it is considered in the lunar module of the Apollo project to 
reduce the temperature extreme of a component (on an antenna) . Another part 
of the heat of fusion experiment is to measure and observe the phase change of 
materials imder a zero-g environment. This experiment will be separated and 
carried under a separate experiment number as a manned flight experiment, 
whereas the radiator portion could be flown as unmanned or manned experi- 
ment. Figure 101 shows the result, of a study which included the possible 
incorporation of the above phase-change experiment in the apollo command 
module. Figure 102 shows a photo of the phase change observation unit of this 
experiment which has been built under a contract to Northrop. In connection 
with this experiment and in support of it, research is carried out in our 
laboratories. The experimental studies include crystal growth, parameter 
studies of various modes of convection, temperature gradients, drop tests 
to simulate zero gravity, differential calorimeter measurements, interface 
between solid and liquid. The emphasis in these experiments is to arrive at an 
understanding of the basic mechanisms involved in crystal growing and the 
solid-liquid interface under one-g and zero-g conditions. 
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FIGURE 98. PORTABLE SOLAR REFLECTOMETER 


188 



HEAT OF FUSION RADIATOR 

ADIABATIC TEST MODEL 


I IB cm 



•“FUSEABLE MATERIAL t-i.jcm 

FIGURE 99. FLIGHT EXPERIMENT FOR MEASURING HEAT OF FUSION 
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FIGURE 100. LAB TESTS OF TEMPERATURE - TIME CURVES FOR 
DIFFERENT RATES OF HEATING PER UNIT AREA 
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FIGURE 101. AAP FLIGHT EXPERIMENT T-005 - APOLLO 

CREW COMPARTMENT 







FIGURE 102. AAP FLIGHT EXPERIMENT T -0 05 - OBSERVATION UNIT 
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ABSTRACTS OF PROGRESS REPORTS 


BUOYANCY EFFECTS ON THERMAL BURNOUT IN 
VERTICAL channels! 


By 

W. R. GambiU 

Oak Ridge National Laboratory 


Abstract 


The effect of buoyant forces on burnout heat flux, which can be signi- 
ficant in low-velocity downflow systems, has been generally neglected. We 
are investigatii^ the subject phenomenon in three related areas: ( 1) burn- 
out of vertical channels blocked at the base and open at the top to a liquid 
Supply, (2) burnout of vertical channels cooled by a low-velocity downflow 
of water, and (3) flow reversal or inversion phenomena in the HFIR {Higfi 
Flux Isotope Reactor) geometry. 

In the blocked-base case, the hypothesis that the critical heat flux 
corresponds to that flux producing sufficient vapor to cause a transition from 
countercurrent Liquid-vapor flow to cocurrent upward flow (flooding) has been 
successful. For a series of 56 tests conducted with round and square 
electrically heated tubes at pressures near atmospheric and with upper pool 
subcoolings of 0 to 142® F, the quantitative expression of this hypothesis may 
be written as follows: 

in which <p is the critical heat flTix, X is the latent heat of vaporization, 
c 

is the equivalent diameter of the channel, and p^ are liquid and vapor 


1. Research sponsored by the U. S. Atomic Energy Commission under 
contract with the Union Carbide Corporation. 
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densities, respectively; g is the gravitational acceleration; and L is the 
heated length. 

Low-velocity downflow tests conducted at atmospheric pressure with 
water have shown that the difference between upflow and downflow burnout 
fluxes Increases monotonic ally with decreasing coolant velocity for velocities 
less than <^1 meter/sec, and that upflow fluxes may be as much as ~7-fold 
larger than those observed with downflow under conditions otherwise the same. 
It has been postulated, but not yet demonstrated, that the minimum critical 
heat flux in a buoyant system will occur with a net liquid downflow ]ust equal to 
the gravitational rise velocity of bubbles in the channel, in which case the net 
body force on such bubbles is zero. 

If the downflow rate in a heated channel is reduced sufficiently, a flow 
reversal or inversion will ultimately occur if the geometry provides a parallel 
flow path. In tests of the HFIR geometry, it was found that the downflow 
velocity just preceding reversal to free upflow could be characterized by the 
ratio of channel-average Grashof to von KSrmSn moduli ( a grouping also 
useful in correlating the results of the low-velocity-downflow burnout experi- 
ments) . 


ULTRAVIOLET RADIATION DAMAGE 


By 

M. L. McDaniel 

Brown Engineering Company, Inc. 


Abstract 


The effectiveness of thermal environment design, as applied to thermal 
coatings presently employed on this nation's space vehicles, depends upon the 
reliability of the optical constants of the material comprising this coating. 

Zinc oxide, exhibiting the most useful solar absorptance to emittance ratio, is 
presently employed as the thermal coating pigment However, the experience 
with this material has been that under the conditions of solar irradiation in 
vacuum, the solar absorptance coefficient is not constant, faulting the thermal 
design. In the present work, the degrading conditions of radiation in space 
have been simulated in an effort to determine the mechanism whereby the 
thermal-optical qualities of this material are compromised. 


194 



PROBLEM AREAS IN THE COMPOSITIONAL ANALYSIS OF 
THE LUNAR SURFACE BASED ON MEASUREMENTS OF 
FAR-INFRARED SPECTRAL EMISS ION OBTAINED BY 

REMOTE SENSING 


By 

J. M Scarborough 
Brown Engineering Company, Inc 


Abstract 


Problems confronting the analyst of far-infrared radiant emission 
from the lunar surface are enumerated and discussed Two simple mathematical 
models incorporatii^ salient features of the radiating process are described. 

One model simulates the radiating properties of a surface composed of 
monolithic grains much larger than the wavelengths of the radiation. The other 
simulates a surface composed of very small particles agglomerated into larger 
clumps which are arranged in an open structure of extreme roughness on a 
scale large compared to the wavelength. 


THE EFFECT OF CRATERS ON LUNAR 
THERMAL RADIATION 


By 

J. P Doty 

Brown Engineering Company, Inc . 


Abstract 


A mathematical model of the lunar surface has been constructed to 
explain the directionality of the lunar thermal radiation. For this purpose a 
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surface with heinlspherieal craters has been assumed, A perturbation on the 
radiation intensity, caused by these craters, is added to Lambert's radiation 
cosine law. The resulting equation is shown to agree with experimental data 
very well. Through the use of this equation, the local' crater density is 
calculated from a general crater distribution of N = NqV where r is the radius 
of the crater 


THERMALLY STRATIFIED SHEAR FLOWS 


By 

Hsing Chuang 
University of Louisville 


Abstract 


The mean velocity distribution in the inner or wall region outside of the 
viscous sublayer of a shear flow is always assumed to be logarithimc if the flow 
field is neutral, but it is approximately expressed as a log-plus-linear profile 
if the flow field is thermally stratified. The extent of deviation from the 
logarithmic profile depends largely on the temperature gradient across the 
region. This is quite understandable from the turbulent momentum transfer 
mechanism. 

In heat transfer computation, it is still a practice to use a logarithmic 
velocity distribution in a neutral flow for the deduction of the turbulent Prandtl 
number. However, both the laboratory and field data indicate that the 
coefficient of the logarithmic term, u.yk, is a fiinction of a thermal stability 
parameter, Richardson number. 
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FREE CONVECTION HEAT TRANSFER TO AIR BETWEEN 
NARROW VERTICAL CHANNELS' 


By 

K, G. Soderstrom and Carlos Reoyo-Sanchez 
University of Puerto Rico 

Abstract 


An experimental model of adjacent fuel plates of the Puerto Rico Nuclear 
Center Reactor was built using the same size plates and spacing as those in the 
Puerto Rico Nuclear Center Reactor. Data have been taken of several variables 
to determine the convection coefficient between the plates and the air between 
the plates in free convection. The results will be used to calculate the 
temperature rise of the fuel plates of the reactor in the case of a sudden loss 
of water under reactor operation at various power levels . 


THE LAMINAR BOUNDARY LAYER WITH HYDROGEN INJECTION 
INCLUDING DIFFUSION DUE TO THERMAL AND 
CONCENTRATION GRADIENTS 


By 

W“S. Lanier and F. Shahrokhi 
The University of Tennessee Space Institute 


Abstract 


When analytical solutions to the problem of hydrogen injection into a 
laminar boundary layer are sought, many of the simplifications which 
characterize boundary layer investigations are no longer valid In the present 
work, hydrogen is injected throu^ a porous wall into a boundary layer 
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composed of nitrogen and oxygen. It is assumed that equilibrium chemistry 
as idealized by the flame sheet model will prevail Further it is assumed 
that a single coihbustion product, water, results from the combustion of 
hydrogen and ojqrgen. The nitrogen present is taken to be an inert species . 
Because of the molecular weights involved, multi-component diffusion must be 
employed instead of the usual Tick's law. Also, there will be an extreme 
temperature gradient which disallows dropping the thermal diffusion term from 
the diffusion velocity equation. In most investigations this term is dropped at 
the expense of accuracy and generality due to the difficulty involved with 
evaluating the thermal diffusion coefficient. 

By assuming a rigid sphere molecular model, the thermal diffusion 
coefficients have been evaluated through the use of Somne expansion 
coefficients. The results for the diffusion coefficients have been specialized 
to the problem at hand but may be easily modified to yield results for an "n" 
component flow field. 


TWO-PHASE FLOW DURING THE FORCED 
CONVECTION ELECTROLYSIS OEWATER^ 


By 

J E. Funk and J. F. Thorpe 
University of Kentucky 


Abstract 


Recent interest in the water electrolysis process for life support 
systems and commercial plants for producing hydrogen and oxygen has created 
a need for information about the two -phase flow which occurs during water 
electrolysis 

This progress report will describe an experimental apparatus which 
will be used to determine the two -phase pressure drop and void fraction 


1. This research is being supported by the National Science Foimdation under 
Grant NSF GK-775. 



distribution during the forced convection electrolysis of an aqueous potassium 
hydroxide electrolyte using stainless steel electrodes at atmospheric pressure. 
The data will be compared with the homogeneous and Martinelli -Nelson models 
of two -phase flow. 

The design and construction of the first electrolysis cell and recirculation 
loop has been completed. Preliminary data have been obtained and will be 
presented, along with a brief description of the planned experimental program. 


D I R ECT I ONAL D I MEN S 1 0 N AL A N AL YS I S 


By 

Dr. R, L Young 

Abstract 


Directional (or vector) dimensional analysis is shown to yield more 
information than that obtained from the conventional dimensional analysis 
Applications of this method to fluid mechanics and heat transfer problems are 
described and discussed. 


DETERMINATION OF THE MOON'S TEMPERATURE AND 
THERMAL PROPERTIES FROM MICROWAVE STUDIES 


By 

J. D. Eggert 

Brown Engineering Company, Inc 


Abstract 


The methods of calibrating microwave receivers for lunar studies 
are described The most consistent and significant experimental results are 
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presented and the importance of this data in drawing conclusions about the lunar 
thermal environment is discussed. Relating the microwave measurements to 
the depth dependence of the Moon's temperature is a central problem in this 
study. The possibility of anal 3 rtioally solving this problem is considered. 


THE EFFECT OF VI BRATIOIM OIM CONVECTIVE HEAT TRANSFER 

IN ENCLOSURES 


By 

E. H. Bishop, C J. Bell, and C, T Carley 
Mississippi State University 


Abstract 


This research task includes both analytical and experimental efforts 
directed at determinii^ the effect of vibration on convective heat transfer in 
enclosures. Published results of the effect of vibration on convective heat 
transfer are apparently limited to cylinders, planes, etc., submerged in an 
infinite fluid. No previous work has been done for a fluid confined in an 
enclosure . 

A mathematical model has been evolved which yields a computer 
solution for velocity and temperature fields. An experimental test cell has been 
constructed for rectangular enclosures and is currently undergoing "shade- 
down" tests on a 2000-pound-force, university-owned, electromagnetic vibration 
system. 

A third portion of the total research effort is directed toward 
determining the effect of intense vibration on fluid properties, especially m 
the cryogenic region. 

The objective of this research is to determine the effect of vibration on 
heat transfer in various geometrical enclosures and to develop a predictive 
ability, for design purposes, for any effects found 
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REPORT ON STUDIES OF TRANSIENTS IN COMPOSITE SYSTEMS 

WITH SURFACES OF CONTACT 


By 

Harold A. Blum 
SMU Institute of Technology 


Abstract 


The major concern of investigators who have been dealing with heat 
transfer across surfaces in contact has been with the determination and 
prediction of the contact coefficient of heat transfer. Most of these studies 
have been carried out under steady state conditions. 

We have been concerned with the heat transfer across surfaces in 
contact in composite systems but under transient conditions . In one situation 
we started with a one-dimensional system consisting of a sink, a source, and 
two metallic cylinders in contact. The system, initially at a uniform 
temperature, was suddenly exposed at one end to the source temperature. In 
this case, the contact conductance coefficient was essentially constant In 
another situation starting with the system just described but initially with a 
steady state temperature distribution, the contact conductance coefficient was 
suddenly increased by a rapid charge of the contact pressure In a third case 
the contact conductance coefficient was changed by suddenly allowing air to 
enter the space which was previously evacuated. 

Both numerical and analytical solutions of the applicable differential 
equations for the three cases have been verified by experiments so that given 
the properties of the metals and the contact conductance coefficients, the 
temperature -time response to either changes in environmental temperature or 
contact conductance coefficient changes can be predicted accurately enough for 
most engineering purposes . The complete description of this study can be 
found in the dissertation of C J. Moore, Jr. (SMU, April 1967) . Some 
aspects of thermal inertia which were predicted by the theory also showed up 
experimentally. 
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The situation in -which one or naore heat sources will be placed on a 
common sink with surfaces in contact has also been considered in our study. 

In this case one would be concerned with the temperature -time response of 
such systems and with the area of thermal influence {that is, the area beyond 
the source in which thermal effects of the source are sensed by the sink) . I'or 
this purpose a numerical solution was developed (two dimensions) for a block 
placed upon a plate. This solution has been checked by ( 1) developing a steady 
state solution and observing the match between the transient and steady state 
solutions as the time becomes large and ( 2) by adjusting the geometry of the 
system for the two dimensional solution in such a way that we could compare 
it with the one dimensional ease. This has been done successfully and some 
runs (computer) have been obtained. 

During the course of this study a passive thermal control device was 
conceived and constructed which made use of a changing contact surface as 
the means for heat transfer control at a given temperature level. 

This study has been supported by a grant from the National Aeronautics 
and Space Administration (NsG71i/44-007-004) . 

EXPERIMENTAL DETERMINATION OF THE CONSTANT 
IN THE CORRELATION FOR NATURAL CONVECTION 
FROM A HOR IZONTAL FLAT PLATE 


By 

Walter Frost and Kenneth L Anderson 
The University of Tennessee Space Institute 


Abstract 


An investigation has been initiated to ascertain the correct value of the 
constant, C, which appears in the correlation for heat transfer by turbulent 
natural convection to a liquid from a horizontal heated flat plate. The standard 
equation given in heat transfer texts [1, 2, 3, 4] is 


Nu 




1/8 
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where C is assigned a value between 0. 13 and 0 17 Recent pool boiling 
heat transfer experiments , [5, 6, V] however, show the value of the constant 
as being on the order of 0.30; a factor of 2.5 difference A closer examination 
of the true value of C is thus needed. 

Preliminary data taken from a 5 x 4| in electrically resistance heated 
plate contained m a 6 x 12 x 12 in glass container gives 0.28 for C, Further 
experiments with different sizes plates and with different fluids are planned. 
Heater plate size variation will discern the length-width ratio necessary for 
true Simulation of an infinite heating surface. The use of dissimilar fluids will 
establish if physical properties of the fluid effect the constant. It is hoped that 
this study will assist in illuminating the reasons underlying the reported 
disagreement in C. 

1. Hsu, S T.: Engineering Heat Transfer. D, Van Nostrand Company, Inc., 
New Jersey, 1963. 

2. Eckert, E. R. G. and Drake, R. M., Jr.; Heat and Mass Transfer. 
McGraw-Hill Company, Inc., New York, 1959. 

3. Jalcob, M : Heat Transfer. Vol. 1, John Wiley and Sons, Inc , London, 
1949 

4 Kreith, F. : Principles of Heat Transfer. International Textbook Co. , 
Pennsylvania, 1958 . 

5. Zuber, N. : Nucleate Boiling: The Region of Isolated Bubbles and the 
Similarity with Natural Convection. International Journal of Heat and Mass 
Transfer. Vol. $, 1963, p. 53. 

6. Bushell, A. L.' Comments on Nucleate Boiling Heat Transfer Correlations 
and Suggestions for Evaluating the Number of Boiling Sites. Royal Air- 
craft EstablishmentrEngland, Technical Note No. Mechanical Engineer- 
ing 391, 1963. 

7. Calculated by the present authors from the data given in; - 
Westwater, J. W. and Kirby, D B., AIChE paper 14, Sixth National 
Heat Transfer Conference, Boston, August, 1963. 
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BI-DIRECTIONAL THERMAL AND SPECTRAL INVESTIGATIONS 
OF POSTULATED LUNAR MATERIALS 


By 

C. D. Cochran 

George C. Marshall Space Flight Center 


Abstract 


The thermal properties of the lunar surface have been measured from 
earth with infrared instruments using an astronomical telescope. From these 
studies, the lunar surface was found to have directional properties and to be 
nonspectral in the infrared portion of the spectrum. This progress report 
presents the results of laboratory experiments with postulated lunar materials 
to study the effect on bidirectional properties of reflectance and emittance, 
particle sizes and mixing of rocks and minerals. 


UNSTABLE NATURAL CONVECTION FLOW IN HORIZONTAL 

CYLINDRICAL ANNULI 


By 

E. H. Bishop, C, T. Carley, and R. E. Powe 
Mississippi State University 

Abstract 


This progress report describes some of the results of an experimental 
investigation of the natural convection of air in horizontal cylindrical annuli 
Specifically, it deals with an unstable oscillatory flow condition recently 
observed by the authors , 
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The flow oscillation is described in terms of dimensionless parameters 
containing the amplitude, period, and wavelength of the oscillations. Empirical 
eq.uations are given which predict the variation of these parameters with annulus 
pressure, Grashof number, and diameter ratio. An equation which gives the 
incipient Grashof number is also presented. 

Work is continuir^ in this area In an attempt to categorize completely 
natural convection flows in cylindrical annuli in terms of the various 
dimensionless parameters pertinent to flows of this type. The present effort 
is directed toward determining the role of the radius of curvature of the 
cylinders on the flow characteristics. 

STUDIES OF SOLAR RADIATION 


By 

Glenn E Daniels 

George C Marshall Space Eli^t Center 


Abstract 


"The Terrestrial Environment Branch, Aero-Astrodynamics Laboratory 
has a special interest in solar radiation at the earth*s surface. This is to 
define solar radiation such that it can be used to determine thermal loads on 
the space v'ehicle on the launch pad and on the components at times when exposed 
to the sun in fabrication, testing, and transportation 

In the document, "Terrestrial Environment { Climatic) Criteria 
Guidelines for Use in Space Vehicle Development," the spectral distribution of 
solar radiation is presented for an extreme value of 1142 watts per square 
meter. The total spectral values versus time of day are also given for the 
extreme and the 95 percentile probability level. These data were based on a 
detail analysis of data from two solar radiation stations; Apalachicola, Florida, 
and Santa Monica, California. 

Future studies are planned to evaluate the diffuse sky radiation for 
extreme values obtained on clear days and for overcast conditions. 
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ENERGY TRANSFER THROUGH EVACUATED 
NONMETAL 1C MATERIALS 


By 

D. W. Tarbell 

Brown Engineering Company, Inc . 


Abstract 


An equation for the apparent theimal conductivity of a quasi -transparent 
medium due to the transfer of energy by radiation has been applied to the case 
of a medium composed of tiny glass spheres in a highly evacuated environment. 
A tentative conclusion is that for spheres with a maximum diameter of 800 
microns, this contribution should be important at temperatures higher than 
200‘’K. 
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APPENDIX A 

ORGANIZING COMMIHEE 


H. L. Atkins 

George C. Marshall Space Flight Center 


Herbert W. Hoffman 
Oak Ridge National Laboratory 


J. Edward Sunderland 
North Carolina State University 


Gerald E. Tanger 
Auburn University 


Reginald I. Vachon 
Auburn University 


Robert L. Your^ 

University of Tennessee Space Institute 
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9:00 AM 

9:15 AM 
10-30 AM 

11:45 AM 

1:00 PM 
1:15 PM 
1:30 PM 

1:45 PM 

2:00 PM 
2:15 PM 


APPENDIX B 
PROGRAM 


Monday, May 8, 1967 


Welcome: R. W. Cook, R&D Operations, MSPC 
SPECIAL LECTURES 


R. B. Bird, University of Wisconsm, "Rheology" 

R, E. Sonntag, University of Michigan 

"Low Temperature Thermodynamics — Hydrogen-Helium 
Systems" 


Lunch 


PROGRESS REPORTS 


W. R. Gambill, Oak Ridge National Laboratory 

"Buoyancy Effects on Thermal Burnout in Vertical Channels" 

L. McDaniel, Brown Engineering 

"Ultraviolet Radiation Damage" 

H. W. Hoffman, Oak Ridge National Laboratory 

"Incipient and Stable Boiling Superheat for Potassium 
Boiling on a Surface Containing Reentrant Cavities" 


J, M. Scarborough, Brown Engmeering 

"Problem Areas in the Compositional Analysis of the 
Lunar Surface Based on Measurements of Par -Infrared 
Spectral Emission Obtained by Remote Sensing" 

J. P. Doty, Brown Engineering 

"The Effect of Craters on Lunar Thermal Radiation" 

H . Chuang, University of Louisville 

"Thermally Stratified Sheer Flows" 
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2:30 PM - Coffee 


2:45 PM - K, G. Soderstrom, University of Puerto Rico 

"Free Convection Heat Transfer To Air Between 
Narrow Vertical Channels" 

3:00 PM - W, S, Lanier and P. Shahrokhi, University of Tennessee Space 
Institute 

"The Laminar Boundary Layer with Hydrogen Injection 
Including Dif&ision Due to Thermal and Concentration 
Gradients" 

3:15 PM - J. E, Funic, University of Kentucky 

"Two-Phase Flow During the Forced Convection Electrolysis 
of Water" 

3:30 PM - R. L. Young, University of Tennessee Space Institute 
"Directional Dimensional Analysis" 

3:45 PM - L, W, Crawford, Tennessee Technological University 

"Fluid Dynamics of a Plane Turbulent Subsonic Air Jet 
Impinging Normally on a Wall" 

4:00 PM - J, Eggert, Brown Engineering 

"Determination of the Moon's Temperature and Thermal 
Properties from Microwave Studies" 

4:15 PM - H. W Butler, West Virginia University 

"Variational Methods in Non-equiiibrium Thermodynamics" 

4:30 PM - E, H. Bishop, Mississippi State University 

"Unstable Natural Convection Heat Transfer in Enclosures" 

4:45 PM - Adjourn 

6:30 PM - Reception 

7:30 PM - Dinner 
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Tuesday, May 9, 1967 


SPECIAL LECTURES 

8 ’30 AM - R, W. Truitt, North Carolina State University 
"Some Aspects of Aerodynamic Heating" 

9:4t5 am - G. B. Heller, Marshall Space Plight Center 
"Thermophysics Research" 

11-00 AM - Lunch 

PROGRESS REPORTS 

12.15 PM - H. A. Blum, Southern Methodist University 

"Report on Studies of Transients in Composite Systems 
With Surfaces of Contact" 

12:30 PM - W, ProstandK. L. Anderson, Umversily of Tennessee Space 
Institute 

"Experimental Determination of the Constant in the 
Correlation for Natural Convection from a Horizontal 
Flat Plate" 

12:45 PM - C. D. Cochran, Marshall Space Plight Center 

"Bidirectional Thermal and Spectral Investigations of 
Postulated Lunar Materials" 

1-00 PM - D. Irey, University of Florida 

"Errors in One Dimensional Fm Solution" 

1-15 PM ^ C. J. Bell, Mississippi State University 

"The Effects of Vibrations on Natural Convection Heat 
Transfer in Enclosures" 

1:30 PM - G. E Daniels, Marshall Space Flight Center 
"Studies of Solar Radiation" 

1:45 PM - D. W. Tarbell, Brown Engineering 

"Energy Transfer Through Evacuated Nonmetallic 
Materials" 
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2:00 PM 


Coffee 


2:15 PM 

2:15 PM 

2: 15 PM 
3:00 PM 

4:00 PM 


GROUP DISCUSSIONS 

I Heat Transfer in Cryogenic Systems - Leader, W. R. 
Gambill, OEJSTL 

n Radiation Heat Transfer in Space - Leader, W. C. Snoddy, 
MSFC 

III Topic to be elected 

General Session: Discussion of Bylaws and Plans for Future 
Seminars 

Adjourn 
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APPENDIX C 
BY-LAWS OF THE 

SOUTHEASTERN SEMINAR OF THERMAL SCIENCES 


article 1. Name 


This organization shall be known as the SOUTHEASTERN SEMINAR ON 
THERMAL SCIENCES, hereinafter referred to as the "Seminar. " The Seminar, 
following the adoption of these By-Laws, shall be denoted as the Fourth 
Southeastern Seminar on Thermal Sciences, and subsequent conferences shall 
be numbered in sequential order. Three seminars having been conducted on a 
"de facto" basis establishes the aforementioned numbering. 


article 2 . Purpose 

The purpose of the Seminar is to stimulate interest in Thermal Sciences 
in the Southeastern area by providing an outlet for technical paper on the results 
of scientific and engineering research of both applied and theoretical nature and 
by encouraging informal exchange of ideas in the field on Thermal Science. The 
Seminar will embrace heat transfer, mass transfer, thermodynamics, energy 
conversion, and fluid mechanics. 

To this end it will arrange for technical sessions of the Seminar to be 
held at approximately two-year intervals or until a change in interval is voted 
by the Policy Committee of the Seminar . The sessions will be held on even 
numbered years unless a change is voted by the Policy Committee It will 
also arrange for the publication of the proceedings of the technical sessions and 
possibly of other scientific research reports. 


ARTICLE 3. Membership 


All interested colleges, universities, and industrial organizations in 
the Southeastern Area of the United States may be members of the Conference. 
The Southeastern Area of the Umted States as used here mcludes Louisiana, 
Arkansas, Kentucl^, West Virgima, Virginia and all states to the east and south 
thereof and the Commonwealth of Puerto Rico. An organization desiring 
membership in the Seminar shall so state, in writing, to the Executive 
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Chairman. Upon receipt of a letter of recognition from the Executive Committee, 
the petitioner becomes a member when: (a) a designated representative from 
the organization is appomted to the Policy Committee, and (b) the Executive 
Committee is notified of the appointment. 


ARTICLE 4. Executive Committee 

A. Composition . The Executive Committee is an appointed body and 
consists of three or four members. There are an executive chairman, an 
executive vice-chairman, and a secretary-treasurer 

B . Duties , The executive chairman has full responsibility for the 
Seminar The Executive Committee is the overall governing body for the 
current technical session of the Seminar and for conducting the business of the 
Conference in the interim between seminar sessions. The Executive Chairman 
is the chairman of the Policy Committee. The secretary-treasurer is 
responsible for collecting and organizing records of executive committee 
actions and correspondence, for passing them on to his successor, and for 
occasionally making permanent disposition of old correspondence, and shall 

be accountable and^responsible for all collected funds, 

C. Appointment of Executive Committee Members. The executive 
vice-chairman is appointed by the host organization for the next succeeding 
technical session and becomes executive chairman for that session. The 
remaimng members of the Executive Committee are appointed by the host 
orgainzation for the current technical session of the Seminar. The current 
technical session is defined as the session for which the newly appointed 
officers are responsible. 

D . Assumption ^ Office, The newly appointed Executive Committee 
takes office at the conclusion of the Business Meeting, as defined m Article 
6D, following the Seminar Session. 


ARTICLE 5. Policy Committee 

A. Composition . The Policy Committee consists of one designated 
representative from each member organization in the Southeastern area and 
the executive chairman of the Seminar. The remaining members of the 
Executive Committee and the immediate past Executive Committee shall be 
ex-officio members of the Policy Committee 
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B . Duties . The Policy Coramittee is responsible for the organization 
and locations of all Semmar Sessions. The members are responsible for 
soliciting papers from their respective organizations. Amendments to or 
changes in the by-laws can be made only by the Policy Committee. The 
secretary -treasurer shall keep a continuing record of its actions and 
correspondence. 

At least three months prior to each Business Meeting the Executive 
Chairman shall request, in writing, the name of the representative to the 
Policy Committee from each member orgamzation. 


ABTICLE 6. Conference Sessions 

A. Selection of Sites for Seminar Sessions . Before a Seminar session, 
an institution desiring to be host to a Conference session should commumcate 
with the Executive Committee. The Executive Committee may also solicit 
invitations from prospective host institutions. The Executive Chairman shall 
inform the Policy Committee of the Seminar of all invitations received, and 
may recommend acceptance. The decision to accept is, however, made by 
majority vote of Policy Committee members present at the Business Meeting. 

If no invitations have been received, prior to the Business Meeting, for 
sponsoring the session following the next session, the new Executive Committee 
shall solicit invitations and a special Policy Committee meeting shall be called 
to make the decision as to where that seminar session will be held. If, at the 
Business Meeting, no site is available for the next session, the Seminar is 
automatically dissolved. 

B. Responsibility for Planning Next Seminar . The actual planning of the 
next Seminar session will be done by officials or faculty members at the host 
institution, which is also finan cially responsible for aU arrangements. The 
host institution is responsible for the formation of an editorial committee, 
which may include persons from other institutions, to select contributed papers 
and invited papers for the Seminar sessions. 

The Policy Committee may provide financial assistance to the next host 
institution, either by an outright grant or by a guarantee of a certain number of 
paid registrants or both, if the Seminar has funds available for this purpose. 

C Publication of the Proceedings of the Technical Sessions of the 
Seminar. The host orgamzation has full responsibility for publishing the 
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proceedings of the technical sessions for which it is responsible. The 
purpose of publication and the guiding principle in the choice of publication 
arrangements is the widest possible circulation of the papers presented at the 
Seminar session. 

D. Business Meeting of the Seminar . In conjunction with each regular 
session of the Seminar, a Business Meeting shall be held. The attendance at 
this meeting shall be limited to members of the Policy Committee, and non- 
participating guests as approved by the Chairman. The agenda shall include: 

1. Report on business transacted since the last Business Meeting, 
including a Financial Statement. These may be duplicated and 
distributed in advance of the Business Meeting and need not be read 
in full at the meeting. 

2. Selection of sites for fixture Seminar Sessions. 

3. Announcement of secretary and executive vice-chairman of 
Executive Committee by the appropriate organizations. 

4. Other business presented by Policy Committee members. 

An opportunity for Committee members to introduce new business shall 
also be provided. 

If procedural questions arise, Roberts Rules of Order shall apply, 
unless the rules are changed by a majority vote of the Committee members 
present at the Business Meeting. 

Interim business may be transacted according to the second procedure 
for amendments as given by Article 7 of these by-laws. 


ARTICLE 7. Amendments 

Two alternative procedures are available for amending these by-laws: 

First Procedure: Proposed amendments to these by-laws shall be made 
by Policy Committee members only, and shall be submitted in writing to the 
Executive Chairman at least three days before the beginning of liie session, and 
copies distributed on the first day of the session to all Policy Committee 
members present. Additional copies of the proposed amendments shall be 
available at the Business Meeting. 
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Amendments so proposed maybe adopted by an affirmative vote of two- 
thirds of the Policy Committee members present and votii^ at the Business 
Meeting. 

Second Procedure : The Executive Committee may propose amendments 
between sessions by mail ballot sent to all Policy Committee members. At 
least thirty days shall be allowed for replies to reach the Executive Committee, 
and at least 50% of the ballots must be returned. Under these conditions the 
amendment may be adopted by an affirmative vote of two-thirds of the ballots 
returned 

Approved at Third Annual SESTS 
Huntsville, Alabama 
May 9, i967 
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Huntsville, Alabama 
May 8-9, 1967 
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University of Tennessee Space 
Institute 

Tullahoma, Tennessee 37388 

Arendale,' W. F. 
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P. O. Box 1247 
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Atkins, Harry L. 

Space Sciences Lab, R-RP-T 
NASA, Marshall Space Flight Center 
Huntsville, Alabama 35812 

Bell, Charles J. 

Mississippi State University 
State College, Mississippi 39762 

Bird, Byron R. 

University of Wisconsin 
Madison, Wisconsin 53706 

Bishop, Eugene H 
Mississippi State University 
State College, Mississippi 39762 

Blum, Harold A. 

Southern Methodist University 
Dallas, Texas 75222 


Caiola, N. J. 

IBM Corporation 
P. O. Box 1250 
Huntsville, Alabama 35807 

Cochran, Charles D. 

Space Sciences Lab, R-RP-T 
NASA, Marshall Space Flight Center 
Huntsville, Alabama 35812 

Chung, Hsing 
University of Louisville 
Louisville, Kentucky 40208 

Crawford, Lloyd 

Tennessee Technological University 
Cookeville, Tennessee 38501 

Daniels, Glenn E 
R-AERO-YT, Bldg 4200 
NASA, Marshall Space Flight Center 
Huntsville, Alabama 358 12 

Doty, J. P. 

Brown Engineering Company 
300 Sparkman Drive 
Huntsville, Alabama 35805 
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North Carolina State University 

Raleigh, North Carolina 27607 
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Eggertj John 

Brown Engineering Company 
300 Sparkman Drive 
Huntsville, Alabama 35805 
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Clems on University 
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Normal, Alabama 35762 

Frost, Walter 

University of Tennessee Space 
Institute 

Tullahoma, Tennessee 37388 

Funk, James E . 
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NASA, Marshall Space Flight Cdnter 
Huntsville, Alabama 35813 

Hoffman, Herbert 
Reactor Division, ORNL 
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Louisiana Polytechnic Institute 
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Klopfenstein, A G. 

IBM Corporation 
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Brown Engineering Company 
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North Carolina State University 
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University of Virginia 
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APPROVAL 

PROCEEDINGS OF THE 
THIRD SOUTHEASTERN SEMINAR 
ON THERMAL SCIENCES 

H. L. Atkins 
R. I, Vachon 

The information in this' report has been reviewed for security- 
classification. Review of any information concerning Department of 
Defense or Atomic Energy Commission programs has been made by 
the MSFC Security Classification Officer. This report, in its entirety, 
has been determined to be unclassified. 

This document has also been reviewed and approved for 
technical accuracy. 



Director, Space Sciences Laboratory 
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